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formed by the Lorentz-noninvariant vacuum 
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Nonlinear electrodynamics, QED included, is considered against the Lorentz-noninvariant external 
field background, treated as an anisotropic medium. Hamiltonian formalism is applied to electro- 
magnetic excitations over the background, and entities of electrodynamics of media, such as field 
inductions and intensities, are made sense of in terms of canonical variables. Both conserved and 
nonconserved generators of space-time translations and rotations are defined on the phase space, 
and their Hamiltonian equations of motion and Dirac bracket relations, different from the Poincare 
algebra, are established. Nonsymmetric, but-in return-gauge-invariant, energy-momentum tensor 
suggests a canonical momentum density other than the Poynting vector. A photon magnetic mo- 
ment is found to govern the evolution of the photon angular momentum. It is determined by the 
antisymmetric part of the energy-momentum tensor. 

PACS numbers: 11.30.Cp, 11.30.Qc, 12.20.-Ds, ll.10.Jj, 13.40.Em, 14.70.Bh. 

Keywords: Lorentz symmetry breaking, Vacuum polarization, Photon Magnetic Moment 



Or 
<D ■ 

(N 
> 

o 

(N 



X 



I. INTRODUCTION 



While relativistic invariance as a symmetry under the 
Lorentz group [SO(3, 1)] is usually an obligatory require- 
ment imposed on a theory, some classes of theories in 
which it is, at least, weakly broken are also of inter- 
est, especially when looking beyond the Standard Model. 
Quantum electrodynamics (QED) in an external classi- 
cal electromagnetic field («^«/) is a clear example of a 
Lorentz-violating theory that may share principal fea- 
tures with other theories of that series. This is a motiva- 
tion for studying it on a very general basis. 

The self-coupling of electromagnetic fields by means 
of the creation and annihilation of virtual charged 
fcrmions makes QED a nonlinear electrodynamics. Like 
any other nonlinear elcctrodynamics-for instance, Born- 
Infcld electrodynamics-QED proposes an interaction be- 
tween a strong classical external field and electromag- 
netic fields that live against its background even when 
these are small. The linearized approximation based 
on the smallness of these perturbations ( "photons" ) will 
be dealt with in this paper. The most important ob- 
ject responsible for the interaction of photons with the 
background in the linearized theory is the vacuum po- 
larization tensor, H^ v {x, x 1 ), calculated in the external 
field. Through this object the gauge sector of QED 
is, in the first instance, provided with a dependence 
on the J^j, structure and therefore on the reference 
frame. Consequently, the photon vacuum seems to be- 



have like an (in general, moving) anisotropic material, 
in which light propagation is strongly modified. Per- 
haps the most remarkable property associated with this 
issue concerns the existence of photon degrees of free- 
dom which arc not in correspondence to the standard 
observable- helicity values-of the respective irreducible 
SO(3, 1)— representations. Instead, the photon propaga- 
tion modes turn out to be closely associated with bire- 
fringent states [1] , and their speeds of propagation differ 
from the speed of light in an empty space-time. Some 
interesting features that occur in the linearized approx- 
imation, when the background is a constant and homo- 
geneous magnetic field, have been predicted. These are 
cyclotron resonance in the vacuum [2] leading to photon 
capture [3], anisotropization, the short-ranging and the 
dimensional reduction of the potential [4] produced by 
a pointlike static charge in a supercritical magnetic field 
\B\ > B c , B c = m 2 c 3 /cfr = 4.42 • 10 13 G, where m and e 
are the electron mass and charge, respectively, 1 and also 
the production of a magnetic field by a static charge (the 
magnctoclectric effect) [5] that takes place in the external 
field, where electric and magnetic fields coexist in paral- 
lel. Beyond the linearized approximation, the important 
effect of photon splitting and merging in a magnetic field 
[6] has attracted much attention. Other consequences of 
the self-interaction of small electromagnetic fields are the 
magnetoelectric effect in QED with an external magnetic 
field [7] and in the nonlinear electrodynamics generated 
by the E/*(l) noncommutative theory [8]. 

Despite the achievements reached in this area, a formal 
treatment of Lorentz Symmetry Breaking (LSB) in non- 
linear electrodynamics has not yet been fully developed. 
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A clear understanding of this theme has paramount im- 
portance in theoretical physics, since the QED vacuum 
in an external field constitutes an ideal laboratory for 
studying the unconventional properties of other Lorentz- 
violating theories encompassed as possible extensions of 
the minimal Standard Model of the fundamental inter- 
actions. Among the candidates appear the Lorentz- 
violating electrodynamics [9-11] and the noncommuta- 
tive field theories [12] with their phenomenologies closely 
related to those described in the preceding paragraph. 
Moreover, since the classification of particles is intimately 
related with the realization of the space-time symmetry, 
the employing of this analysis could lead to new insights 
on plausible phenomena in which other spin representa- 
tions like axions (spin 0) [13] and gravitons (spin 2) [14] 
arc coupled to & . 

Of course, an experimental confirmation of all 
these processes strongly depends on the external field 
strengths, whose current laboratory values are much 
lower than the critical field. As a consequence, all pre- 
dicted phenomena remain elusive and far from being de- 
tectable. Even so, their studies still find a high mo- 
tivation in light of upcoming laser facilities [15, 16] that 
will achieve the unprecedented level of |-E7| ~ 0.01 — 0.1 E c 
with E c = m 2 /c = 1.3 • 10 16 V/cm. In addition, some 
evidence points out the possible existence of ultrahigh 
magnetic \B\ 3> B c and electric \E\ 3> E c fields in the 
surfaces of stellar objects identified as neutron stars [17] 
and strange stars [18], respectively. In such scenarios a 
pronounced LSB is expected and most of the quantum 
processes described above could have significant astro- 
physical and cosmological interest. 

Inspired by the importance associated with LSB, we 
make an attempt to fill some gaps in this topic. Our 
main purpose is to analyze how the vacuum polarization 
effects modify the Lorentz and Poincare generators. The 
results presented in this work are based on the Poincare 
invariance of the photon effective action as a functional 
of the background field. This implies that those Lorentz 
transformations which leave the external field invariant, 
together with the space-time translations that do not af- 
fect the external field either, as long as it is time and 
space independent, provide the residual symmetry sub- 
group of the anisotropic vacuum, while the full Poincare 
group remains the group of broken symmetry. To de- 
fine the space-time translation and rotation generators, 
although only a part of them is conserved, we appeal 
to the context of the Noether theorem and then intro- 
duce them into the framework of the Hamiltonian for- 
malism, which requires us to impose constraints associ- 
ated with the gauge invariance of the theory, and thus the 
Dirac brackets. The general aspects related to the con- 
strained Hamiltonian dynamics were developed by Dirac 
[19] and have been applied to several problems in Quan- 
tum Field Theory [20] , including the analysis of Poincare 
invariance in Yang-Mills theories quantized in noncovari- 
ant gauges [21] (for the Coulomb gauge, sec Rcf. [22]). 
We do not know whether the Hamiltonian formalism was 



ever applied to electrodynamics of an anisotropic medium 
or whether the characteristic entities of the latter are 
known in terms of canonical variables, but in any event 
the present exploitation of this formalism results in some 
interesting features. Among them is the distinction be- 
tween the generating function of infinitesimal canonical 
transformations of spatial translations, that most natu- 
rally turns out to bc parallel to the wave vector in each 
cigenmode, and the Poynting vector that points in the di- 
rection of the energy propagation and the group velocity. 
The Dirac commutation relations between the space-time 
generators leave intact the 50(3)— algebra of the angular 
momentum, although this symmetry group is broken, but 
the rest of the Poincare algebra commutation relations 
that include transformations of the time are distorted. 
The Hamiltonian equation of motion for the nonconserv- 
ing angular momentum indicates that its time evolution 
is determined by the interaction of the external field with 
a certain magnetic moment (that also contributes to the 
Hamiltonian) emerging in the special frame in which a 
pure magnetic field is present. Such a quantity is inter- 
preted as the magnetic moment of the photon; i.e., the 
magnetic moment of the anisotropic medium that is not 
only polarized, but also magnetized by the photon field. 
The possible connection of the photon magnetic moment 
with the same notion introduced previously by Perez Ro- 
jas and one of the present authors (Villalba-Chavez) [23] 
(see also Refs. [24, 25]) will be discussed in a due place 
below. It looks at attribute of the photon interaction 
with an anisotropic medium. 

We organize the manuscript as follows. In Sec. II A, we 
recall some basic features of the photon propagation in an 
external field based on the relativistic covariant formal- 
ism introduced in Ref. [1] that involves diagonalization 
of the polarization tensor and analysis of its eigenvalues 
and eigenvectors. In Sec. II B, we consider the contribu- 
tion of the polarization operator to the effective nonlocal 
action (the generating functional of irreducible vertices 
[26]) to give it the form of the action of the equivalent 
linear anisotropic medium with time and space disper- 
sion. The tensor decompositions and principal values of 
dielectric, e^, and magnetic, permittivities-related 
to the special class of Lorentz frames where there is only 
a magnetic, or only an electric, external field-are pre- 
sented in terms of the polarization operator eigenvalues 
in an approximation-independent way. We point out the 
uniaxial character of the vacuum in these frames. This 
statement holds true also for the most general combi- 
nation of constant electric and magnetic fields in the 
Lorentz frames, where these fields are parallel. In Sec. 
II C, for the same general external field, we obtain the 
covariant decompositions of the global coordinate trans- 
formations that leave the external field intact and thus 
make the invariance subgroup of our problem. 

In Sec. Ill and thenceforth, we confine ourselves to 
the local approximation of the effective action, when it 
does not include field derivatives as the functional argu- 
ments. This limitation will make the subsequent devel- 
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opmcnt of the Hamiltonian formalism straightforward. 
The local approximation corresponds to the small four- 
momentum (infrared) limit of the polarization operator 
and to the frequency- and momentum-independent ten- 
sors Eij and Hij. The latter are expressed in terms of 
derivatives of the effective Lagrangian over the field in- 
variants for the cases reducible to a single field in a spe- 
cial frame (magncticlike and electriclike cases). In Sec. 

III A, based on the Noether theorem, we define a nonsym- 
metrical, but gauge-invariant, energy-momentum tensor 
of small perturbations of the vacuum that satisfies the 
continuity equation with respect to only one of its ten- 
sor indices. The antisymmetric part of this tensor will 
become, in what follows, responsible for the nonconser- 
vation of those Lorentz and spatial rotation generators of 
canonical transformations constructed using the energy- 
momentum tensor, which do not leave the external field 
invariant. We stress the difference between the momen- 
tum flux vector (canonical momentum), whose direction 
is parallel to the wave vector according to Appendix C, 
and the energy flux Poynting vector, whose direction co- 
incides with that of the group velocity and the center- 
of-mass velocity of eigenmodes, as demonstrated in Ap- 
pendix D. The constrained Hamiltonian formalism serv- 
ing the dynamics of small perturbations over the back- 
ground field is presented in Sec. Ill B following the pro- 
cedure well elaborated in gauge theories. It is the electric 
induction of the perturbation, and not the field strength, 
that comes out as a variable, canonically conjugated to 
its three- vector potential. With the Coulomb gauge con- 
dition extended to the problem of the anisotropic vac- 
uum under consideration, the Dirac brackets are defined 
as performing the infinitesimal canonical transformations 
in the phase space. In Sec. IIIC, referring again to the 
Noether transformations, we define the conserved and 
nonconserved components of the angular momentum and 
the Lorentz boost, express them in terms of the canoni- 
cal variables, and find their Dirac brackets with the fields 
and inductions. In Sec. IV, the Hamiltonian equations 
of motion for these quantities are given (Sec. IV A), 
as well as the set of Dirac commutators for the gener- 
ators of space-time translations and rotations, defined 
above, that substitute for the standard relations of the 
the Poincarc algebra in the present case the vacuum in- 
variant not under all Lorentz and space rotations (Sec. 

IV B and Appendices A, B). In Sec. IV C, we dwell on 
the algebra of the space-time invariance subgroup and 
define its conserved Casimir invariants. 

In Sec. V, we deal with the magncticlike external 
field. The magnetic moment 9A. of the photon propagat- 
ing over the magnetized vacuum is analyzed. It appears 
as an entity that governs the evolution of the photon an- 
gular momentum ^ in the magnetic field B following 
the equation of motion d J? /dx° = 2fW x B, and con- 
tributes into the photon energy as — $4. ■ B. A further 
step in our understanding of this quantity is given by 
showing its connection with the optical tensors of our 
problem. In the large— \B\ region, the photon magnetic 



moment treated following the one-loop approximation of 
quantum electrodynamics depends quadratically on the 
photon electric field alone. Its appearance may be un- 
derstood as another manifestation of the magnetoelec- 
tric effect [5, 7] in QED, known also in noncommutative 
electrodynamics [8]. 

In Sec. VI, we write down the coefficient tensor cus- 
tomarily used to serve the gauge sector in the general 
Lorenz-violation approach to a [/(l)— invariant theory, 
as it follows from the general covariant decomposition 
of the polarization tensor in a magnetic field found in 
Ref. [1]. We establish that the coefficient tensor is not 
double traceless, contrary to what is assumed in the 
above approach. The double trace is physically mean- 
ingful as being connected with the magnetic and electric 
permeability of the magnetized vacuum. We express the 
condition of the absence of birefringence in terms of field 
derivatives of the effective Lagrangian. Finally, we esti- 
mate the values of the magnetic field likely to produce the 
Lorenz violation that would be equivalent to the Lorentz 
violations intrinsic in the vacuum and detectable using 
experimental devices of present-day sensitivity. These 
magnetic fields arc too large to make a realistic cosmic 
background. 

We present our concluding summary in Sec. VII, while 
the essential steps of many calculations have been de- 
ferred to the appendixes. 

II. LORENTZ SYMMETRY BREAKING: 
GENERAL ASPECTS 

A. The photon effective action 

In the presence of an external field ^(x) = —\& llv '£ v 
with a constant field strength = d^stfy — d^g/^, the 
action which describes small-amplitude electromagnetic 
waves a^ix) over a constant background field reads 

S = -\J Tf^x + T, (1) 

where = d^a v — <9„a M , and V is called the effective 
action, connected as T = J £,d 4 x with the effective La- 
grangian £. Effective action may be expanded in powers 
of the small field potential as 

T = \j d4xdix ' ati ( x ) u ^( x > x 'W) aU ( x ') + • ■ • (2) 

where + . . . stands for higher-order terms in whereas 
n M1/ (a:, is the second-rank polarization tensor re- 

lated as 

9^{x,x'\^) = [n VfiV - d^d v \ 6^(x' - x) 

+n lu ,(x,at\j*) (3) 

with the inverse photon Green function S$~^(x, x'\srf). 
Hereafter, the metric tensor rf' v has the signature ++H — 
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with rj iL = if" = T] da = —r] = 1; the electric held is 
given by E l = J£" 0i , whereas the magnetic field is defined 
by B l = l/2e ijk 3? jk . The tensor ^ v = l/2e^P a ^ pa 
(with e AIlyptT being the fully antisymmetric unit tensor, 
£ 1230 _ represents the dual of The field in- 

variants are 3 = 1/4J^„J^ = (B 2 - E 2 )/2 and 
= 1/43?^^ = E ■ B = 0. For the special case of 
an external field mainly dealt with in the present work, 
we shall choose © = 0. Then the corresponding external 
field will be referred to as magneticlike (if 3 > 0) or elec- 
triclikc (if J < 0), because in either case a Lorentz frame 
exists wherein the field is purely magnetic or electric, re- 
spectively. 

According to Eq. (2), a photon can interact with the 
external field through the vacuum polarization tensor 
TL iiv {x,x'\j^). In this context, the QED Schwinger- 
Dyson equation for the photon field a fJi (x) is given by 



d i x'@-Z{x,x'\s/)a v {x') = Q. 



(4) 



The external field strength is independent of the space- 
time coordinates; therefore, the polarization tensor, as 
a gauge-invariant quantity, should correspond to a spa- 
tially homogeneous optical medium whose properties do 
not change with time. This is provided by the trans- 
lational invariance of 11^^: it depends only on the co- 
ordinate difference II M! ,(x, x'\srf) = H^(x — x'\g/) [27]. 
In this Fourier transform converts Eq. (4) into a 

linear homogeneous algebraic equation given by 



IVOIO] a v (k) = 



(5) 



with 



iWfcl^H / ^ v {x-x'\,^)e- lk ^- x 'U\x-x'). (6) 



To understand what follows, it is necessary to recall 
some basic results developed in Refs. [1, 27]. In the pres- 
ence of a constant magneticlike or electriclike field, the 
four eigenvectors of the polarization operator b^ are 
known in a final, approximation-independent form. In 



(4) 



A',, 



addition to the photon momentum four-vector b 
(its zeroth component k° being the frequency u>), the 
three other mutually orthogonal four-transverse eigen- 



vectors bp are 



b^_=k 2 ^ x k x -kp(k^ 2 k), 

^yu\k ,(3) _ ^fi\k^ 



" fi •• " fj,\ 

S ;„A 

b {2 - 

11 (k^ 2 ky/ 2 



b (3) =. 
M {-k& 2 ky/ 2 ' 



(7) 



k^b 



(A) 



for A = 1,2,3. (The eigenvectors relating 



to the most general case, J / 0, © ^ 0, are writ- 
ten in Refs. [1, 27] and [28, 29]). We remark that 

bjL fulfills both the orthogonality condition, bp X \^ x '^ = 
gxx (|/A)^ ^ an( j completeness relation, 



km 2 ' 



(8) 



Note that from b^ x \ one obtains the fundamental scalars 



k 2 
k& 2 k 

23 



= Z\ + 2 2 , 

and Z2 = 



k,^ 2 k 
23 ' 



(9) 



The last two scalars acquire simple forms in a special 
reference frame where the external field is purely mag- 
netic (if 3 > 0) or purely electric (when the opposite 
inequality holds). The same equations hold in reference 
frames that are moving parallel to the external field. For 
the magnetic background 3 > 0, one finds that Z2 = k\ 
and z\ = k 2 — lu 2 . On the contrary, if the electric field is 

considered, z 2 = k 2 — uj 2 , whereas Z\ = kj_. The previous 
relations involve the vectors k± and ku , which denote the 
components of k perpendicular to and along the exter- 
nal field, respectively. Henceforth, boldface letters will 
designate the spatial part of our four-vectors. 

Besides the creation of the fundamental scalars, i.e., 
Eq. (9), the vectorial basis b^ A ^ is suitable to express the 
vacuum polarization tensor in a diagonal form 



4 b w bl x) 

IV =^«X(*l,*3,ff) jj^jl 



(10) 



A=0 



Here x\ denotes the eigenvalues of the vacuum polariza- 
tion tensor 



n;bW = x (fc) 



1,2,3,4, 



which define the energy spectrum of the electromagnetic 
waves and poles of the photon propagator. Owing to 
the transversality property (fc^n^ = 0), the eigenvalue 
corresponding to the fourth eigenvector vanishes identi- 
cally [x^ = 0]. Substituting Eq. (10) into Eq. (5) and 
using the orthogonality condition, we find its solutions 
in the form of a superposition of eigenmodes given by 
a f _ l (k) = Fx5(k 2 — x\)b^ X \ where F\ are arbitrary 

functions of k. According to the latter, three nontrivial 
dispersion relations arise: 



k 2 = x\ (z 2 ,2i,3) , A = 1,2, 3, 



whose solutions can be written as 



^ = fcf + A(fci,3). 



(ii) 



(12) 



The term /a (^,3) arises as a sort of dynamical mass. 
Due to the gauge invariance condition xa(0,0,3) = 
there always exist two (out of three) solutions with 
/a (0,3) = that correspond to photons whose rest en- 
ergy is zero, and the number of polarization degrees of 
freedom is two. Massive branches of all the three polar- 
izations, f\ (0,3) 7^ 0, A = 1,2,3, may also exist. For 
more details, we refer the reader to Refs. [5, 28]. 

Moreover, by considering ~ b^ as the electromag- 
netic four-vector describing the eigenmodes, we obtain 
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the corresponding electric and magnetic fields of each 
mode in the special frame, provided that $ > : 

eW ~i(u/ A V A > -fcb°( A >) and ~ ik x b^ . (13) 

Up to a nonessential proportionality factor, they are ex- 
plicitly given by 



~ —in±uj, &W ~ — iA; x n_i_, 



^~-ifc ± fe||/(fejf-« a ) 



2 ,,,2U/2 J2)^_ 



in||(fcf-w 2 ) 1/2 - 



6( 2 ) xn||)/(A; 2 -a; 2 ) 1 /2 ) 



„(3) 



zw(n^ x nil), 



6|| 3) ~ —in\\k±, by ~ m_i_fc||. 

Here, ri|| = fc||/|&||| and n = k±/\k±\ are the unit vec- 
tors associated with the parallel and perpendicular direc- 
tions with respect to the magnetic field B in the special 
frame in which the external electric field E vanishes iden- 
tically. 



B. The vacuum as an anisotropic medium within a 
linear optics approximation 

In classical electrodynamics, the Maxwell action of 
a linear continuous medium with the dielectric tensor 
Eij(k,uj) and the magnetic permeability tensor fiij(k,ui) 
is given by the expression (quadratic in the fields) 

f 1 1 
S = / L = -d(—k, —lj) ■ e(k,ui) - -h(—k, —u) ■ b(k, uj). 

(14) 

Here e, b are connected by means of the relations 



d(k,uj) = e(k,oj) ■ e(k,u), 
h(k 7 uj) = /x -1 (A;,w) • b(k,ui). 



(15) 



where the double-sided arrow denotes a tensorial quan- 
tity. In this context, *e ■ e = e^ej and /2 _1 • b = fi^bj. 

The optical properties of an anisotropic medium de- 
pend primarily on the symmetry of its tensors £y and 
|Ujj . In an uniaxial medium, one of the principal axes of 
Eij and iiij forms the "optical axis" . In what follows, we 
denote the principal values of and /i.y relating to this 
axis as em and /xm, respectively, and the values relating to 
the plane perpendicular to the optical axis as e±_ and fj,± , 
respectively. With this in mind, the Maxwell Lagrangian 
[Eq. (14)] acquires the following form: 



L=Ue^\ 2 -—\b A 
2 I fx± 



■ £ \\\ e \\ 



/'II 



>\ 



(16) 



where we have decomposed e = (e_i_,e|i) and b = (6j_,6||). 
Here the symbols _L and || refer to the optical axis as well. 

Let us consider now the quadratic part of the effective 
action corresponding to the dynamical gauge field sec- 
tor of QED in an external field. Substituting Eq. (10) 
into Eq. (2) and making use of Eq. (8), we find for 



the integrand of S [Eq. (1)] in momentum space 

_ - f» v f + £ 



1 



Sf=~0^f ltv . (17) 
Here the second-rank antisymmetric tensor reads 

= (l - /' 



1 Xi — X3 

2 k^ 2 k 



2 k^ 2 k 



(18) 



The Lagrangian ££ was written in Rcf. [28]; its small- 
momentum form [corresponding to Eq. (47) in the next 
Subsection] is present in an earlier paper [30] . 

The expression above is valid for both magneticlike 
and electriclike cases (J ^ 0, = 0) and defines the cor- 
responding induction vectors according to the following 
rule: 



d* = O kh = d^/de 1 



(19) 



where e(x) = Vao(x) — doa(x) and b = V x a are the av- 
eraged (classical) electric and magnetic fields associated 
with the electromagnetic wave. With these definitions, 
the Maxwell equations have the recognizable form 



Vd = 0, V-6 = 0: 



V xe- -— V x h - 9d 

dx° ' dx° 



(20) 



The above expressions allow us to obtain the most gen- 
eral structure of the dielectric and magnetic permeability 
tensor. To derive it, we first note that our effective La- 
grangian 3f in Eq. (17) acquires the structure of Eq. (14) 
as long as it is expressed in terms of the induction vectors 
d, h and the electric e and magnetic b fields of the small 
electromagnetic waves. Likewise, the optical tensors can 
be defined as in Eq. (15). In fact, by considering the 
following relations, valid in the special frames 



one can express them for the magnetic external field 
e»j(*,w)= 1-TaK" 1 



(21) 



k 2 J 13 k 2 -u 2 B 2 



X! - x 3 BiB, 



(22) 



It is notable that the components of the three-momentum 
vector A; do not take part in forming these tensors; only 
components of B do. This feature is not typical of crystal 
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optics with spatial dispersion and can be attributed to 
the explicit exploitation of the gauge invariance laid in 
Eq. (18). The eigenvalues of matrices in Eq. (22) [the 
principal values of the electric and (inverse) magnetic 
permittivities] are 



_! Xl X X , x x - x 2 



fc 2 fcjj-w 2 ' 



(23) 



-1 1 Xi - X 3 



The values £|| and /^| correspond to the eigenvector di- 
rected along the external magnetic field B, which there- 
fore makes it the direction of the principal optical axis. 
The values £j_ = /i^ 1 correspond to the eigenvectors di- 
rected transverse to the external magnetic field. The 
principal values of Eq. (23) are rotational scalars and 
depend upon direction: their arguments are the fre- 
quency u> and the scalar product B ■ k, combined into 
z\ = bJ 2 — k 2 , z 2 = kj_. We also find a similar result 
for an electriclike background (5 < 0). In this case, the 
optical axis is determined by E and 

(24) 



^ X (*,«)=(l-^) <% + 



Xj - x 2 EjEj 
k 2 E 2 



From these tensors, we obtain that e± and have the 
same structure as in Eq. (23), whereas 



£|| = £j 



Xi - x 3 



k 2 

II 



My 1 = Va 1 



Xi - x 2 



(25) 



The results given in Eqs. (22)- (25) point out that in the 
presence of an external magnetic or electric field, the vac- 
uum behaves like an uniaxial anisotropic material. Note 
that the procedure shown in this section is independent 
of any approximation made in the calculation of the vac- 
uum polarization tensor. However, it is only valid in 
a class of special frames, in which the external field is 
purely magnetic (when 25 = 0,^ > 0) or purely electric 
(when 25 = 0,^ < 0). In a general Lorentz frame, an 
electric (magnetic) component is added to the primar- 
ily purely magnetic (electric) field, as produced by the 
Lorentz boost. Hence, the statement that the vacuum 
is uniaxial is no longer true in that frame, because the 
second axis is specialized by the direction of the added 
componcnt-or, in other words, by the direction of the 
motion of the reference frame with respect to the spe- 
cial frame. Therefore, the vacuum is a biaxial medium 
that can be rendered uniaxial by an appropriate Lorentz 
transformation. (In the case of a material anisotropic 
medium, which is uniaxial in its rest frame, we can also 
state that it becomes biaxial if it moves with respect to an 
observer, the direction of motion specializing additional 
direction in the frame of the observer.) 



The same statements are readily extended to the case 
of a general external field with both the field invariants 
different from zero: 25 ^ 0, $ ^ 0. In this general case, (a 
class of) special frames exist, where the external electric 
and magnetic fields are mutually parallel, their common 
direction specializing the principal optical axis in such 
frames. The point is that the diagonal representation for 
the polarization operator [Eq. (10)] remains valid, the 
only reservation being that now the eigenvectors in 
it are not just the vectors of Eq. (7), but linear com- 
binations of them [29]. A representation analogous to 
Eq. (18) can be written in that case, and the principal 
values substituting for Eq. (23) again depend on the same 
combinations of momenta uj 2 — fcjj and k\ , where now the 
designations || and _L mark the directions parallel and or- 
thogonal to the common direction of the external fields. 



C. The vacuum symmetry subgroup ISOa(3, 1): 
Covariant decomposition of transformations 

Actually, the anisotropic character of the medium, 
equivalent to that of the vacuum with an external field, 
arises due to the Lorentz and rotational symmetry break- 
down [which is not manifest in the Maxwell Lagrangian 
of Eq. (14), because it relates only to the rest frame 
of the medium and does not reflect its spatial symme- 
try]. The Lagrangian [Eq. (17)] is not Lorentz and ro- 
tational invariant, because it contains an external tensor 
responsible for the external field. So, to keep it invariant, 
one should transform the external field together with the 
photon field. Correspondingly, the explicit forms of the 
scalars z 2 and Z\ , as well as Xj , when these are expressed 
through the photon momentum components, depend on 
the reference frame. However, the Lagrangian [Eq. (17)] 
turns out to be invariant under those space-time trans- 
formations which leave the external field intact. Thus, 
bearing in mind the translational invariance of our prob- 
lem, the proper inhomogeneous orthochronous Lorentz 
transformations relating to the symmetry group of an 
anisotropic homogeneous vacuum occupied by an exter- 
nal space- and time-independent classical field must fulfill 
the conditions 



,^=AP \\^ pa , detA = l, A°>0. 



(26) 



The set of pairs {e, A} satisfying Eq. (26) form a subgroup 
of the Poincare group [I SO (3, 1)] which will be referred to 
as the "Amputated Poincare Group," ISOa(3, 1)- Also, 
A 6 SO A (3,l), where SC^ (3,1) is called the "Ampu- 
tated Lorentz Group." Due to Eq. (26), the infinitesimal 
Lorentz transformation associated with our problem can 
be written as [311 



A". 



6". 



with 



(27) 
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where and £' are real infinitesimal parameters. Note 
that the description above is not restricted to a magnet- 
iclike or electriclike field. On the contrary, it holds for 
any other constant and homogeneous external field con- 
figuration with the second field invariant, © ^ 0, and it 
emerges in any other nonlinear electrodynamics different 
from QED. 

We exploit Eq. (27) to obtain explicitly the struc- 
ture of those Lorentz transformations that are associ- 
ated with our problem, i.e., A e 50^(3,1). Before do- 
ing this, we redefine the group parameter in Eq. (27), 
& = ■&/(Af+ + A/1) and £ = f/(A/+ + A/1), where 

1/2 

A/± = [G? 2 + <5 2 ) 1/2 ± $] are eigenvalues of the ex- 
ternal field tensor & . In addition, we express the finite 
Lorentz transformation as an exponential of a matrix ar- 
gument 



A = exp 



M++N- 7V+ + A/l 



(28) 



defined by its series expansion. 

Substantial simplifications can be achieved by the in- 
troduction of the matrix basis 



whose elements fulfill the following properties: 



(29) 
(30) 



Z+Z_ = 0, Z 2 " = z. 



r2n 



(-1) 



n — 1 -2-2 



2 7 

j-271+1 



2)1+1 



= z 4 



(-l) n Z_ 



(31) 



With these details in mind, a finite Lorentz transforma- 
tion belonging to 50^(3, 1) decomposes according to 

A = Z_ sin tp + Z + sinh C, — Z 2 cos tp + Z 2 . cosh (, (32) 

where the the arguments of the trigonometric and hyper- 
bolic functions are 



c = 

<p = 



A/+£ 



A/ltf 



A/+ + A/1 A/V+A/l 
A/l£ A/"+tf 



A/1 + A/1 A/l+A/1' 



(33) 
(34) 



In contrast to the cases considered in Subsections A and 
B, the covariant decomposition of the Lorentz transfor- 
mation [Eq. (32)] deduced here relates to the most gen- 
eral case of a constant and homogeneous external field 
with both its invariants ©, 5" different from zero. (There- 
fore it remains valid in the crossed field system (3=^ = 
and \E\ = \B\). For the magnetized vacuum (0 = 0, 
5 > 0), the variables in Eq. (32) become ( = £ and 
tp = whereas in an electric background (© = 0, 5 < 0) 
they turn out to be Q = and tp = £. 



Now, the explicit structure of A in terms of the exter- 
nal electric E and magnetic B fields is rather complicated 
in a general Lorentz frame. However, it becomes simpler 
if one considers the external field configurations inherent 
to special frames, where the vectors E and B are paral- 
lel, and directed, say, along the z axis. In these special 
frames, the transformations of Eq. (32) have a simple 
sense of rotation about the axis z by the angle tp, and 
the Lorentz boost along this axis is parametrized by the 
angle £. 

Considering the general case 25 =^ 0, $ ^ in the 
special frame where B = (0,0, B z ), E = (0,0, E z ), we 
found that 

=^\{V)@ X AC) = E> 50(2) ®B 50(M) . (35) 

The expression above involves the matrices 

= D 50(2) © 1 and ^(C) = l ©D so(14) , (36) 

where the © denotes the usual direct sum of matrices, and 
II denotes a 2 x 2 identity matrix. Here both B so ( 2 ) and 
j)SO(i,i) are two-dimensional representations of 5*0(2) 
and SO (1, 1), respectively. Explicitly, they read 



D SO(2) = 



cos tp sin tp 
— sin tp cos tp 



D 



SO(l,l) _ 



cosh ( sinh ( 
sinh £ cosh ( 



Both matrices conform two independent Abclian- 
invariant subgroups of 50,4(3,1), which emphasizes its 
nonscmisimplc structure. Moreover, the group parame- 
ter space is the product manifold 5^ x 5ftW, where 5^ 
is the circle and tft^ the real line. Therefore, the topol- 
ogy of 50,4(3, 1) is the surface of a circular cylinder, a 
manifold which is neither compact nor simply connected. 
(Indeed, it is infinitely connected.) 



III. GENERATORS OF SPACE-TIME 
TRANSFORMATIONS 

The constrained Hamiltonian formulation of the elec- 
tromagnetic field in an anisotropic material remains 
poorly developed in the literature; we are not aware of 
previous Hamiltonian treatment for these optical media. 
In this section we shall see how this method can be im- 
plemented in QED with an external constant electromag- 
netic field. To pursue this analysis we first note that the 
effective Lagrangian Jz? in Eq. (17) is not a local func- 
tion, since it depends on momenta in a very cumbersome 
way. The canonical formalism in its standard form is not 
applicable to nonlocal situations. For this reason we will 
restrict ourselves to the infrared approximation, fc M — > 0, 
in which case the effective Lagrangian Jzf becomes a local 
function on the photon field a^x). This approximation 
corresponds to anisotropic media with no spatial or fre- 
quency dispersion. In our case, it becomes actual in the 
region of very strong external fields, where the external 
field dominates over the photon momentum [32] . 
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Consider the static (u> = 0) case. When k± = 0, we 
have from Eq. (23) 



£||(fc||,0) = l 



^2 



When uj = 0, fcii = 0, we have that 



Mjf A (fcx,0) = l-^ 

Owing to the degeneracy property [5], >t\ = x 2 at 
k?, = 0, we also that 



(37) 



ex(fci,0) = A*I 1 (fc±,0) = l- ^ 



j, fcii =0 



So, in the infrared limit fc^ — > 0, the quantities in 
Eq. (23) coincide with the permittivities defined for 
that case in Ref. [28]-according to the correspondences: 

e± e tr , (J>± o Mtr> £|| ^ eiong, M^ 1 /«tr~ and 
responsible for screening charges and stationary currents 
of special configurations. 

In the infrared limit the eigenvalues of the vacuum po- 
larization tensor can be expressed in terms of the first 
and second derivatives of an effective Lagrangian £ [con- 
nected with the generating functional F of irreducible 
many-photon vertices in an external field, as pointed out 
below Eq. (1)] over the constant field with respect to the 
corresponding external field invariants £ and C5 (sec de- 
tails in Ref. [28]): 



(38) 



where % = d£/d$, £«© = <9 2 £/<9© 2 and % s = 
d 2 £/d$ 2 with © set equal to zero after differentiation. 
It follows from Eqs. (38), (22), and (24) that 



Eij - (1 - £ s ) 8ij + ZeeBiBj , 
^ = {1-^)5^-^13^ ' 

eij = (1 - £y) 5^ + £,$d E i E j 



(39) 



For the magncticlike case J > this is equivalent to 



Ml 1 = i - fia 



£ |l =f-£ 5 + 2S£ 0l5 , 



(40) 



The following relations were established in Ref. [28] on 
the basis of causality and unitarity principles, valid for 

1 - £3 > 

1 - £ s + 2$£ ee > 0, £ @& > (41) 
1 - £5 - 2££ 5S > 0, £ S5 >0 

which guarantees the consistency of the theory. 



A. The energy- momentum tensor 

The symmetry reduction by the external magnetic 
field does not alter the translational group embeded in 
ISO(3, 1). Therefore, for a photon, the spacetimc con- 
figuration with an external classical field is translation 
invariant. To find in the local approximation the asso- 
ciated Noether current of the electromagnetic radiation, 
let us first insert Eq. (38) into Eq. (18) to find 



QP.V 



1 



f pX ^ P x 



(42) 



The substitution of this tensor into 5£ [Eq. (17)] de- 
fines the Lagrangian [28, 30] of the small-amplitude, low- 
frequency, long- wave electromagnetic field a^{x): 



= -i(i-%)/^/^ + ^«« I 

+^%3 {r^f 



(43) 



The corresponding conserved stress-energy tensor is ob- 
tained from it by following the Noether theorem in the 
standard way: 



dse 



Here 



d v a x . 



(44) 



(45) 



in accordance with Eqs. (17)-(19). The antisymmetricity 
of this tensor is owing to the gauge invariance manifesting 
in that the Lagrangian contains only the field tensor. 
Substituting Eq. (43) into Eq. (44) we obtain 

x f e «r + j£6 6 (r e A e ) 2 + ^£35 {r^a e f}. 



(46) 

Let us define the tensor 

0/w = c^^v + 0» x r x - fa v . (47) 

related to Eq. (44) as 

T^ v = - d x K x ^ (48) 

where K XflL/ = —0^ x a v is antisymmetric in its first two 
indices, while the electric current 



r = 9^ 



d{du,a u ) 



(49) 
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disappears on equations of motion = 0. It conserves, 
d^j^ = due to antisymmetricity of the tensor 0^ x 
[Eq. (45)]. When taken on the equations of motion, the 
tensors T^ v and 8'"' coincide up to a full derivative. 
However, the latter depends only on the field strengths 
of the small electromagnetic field f a p, and not on its 
four- vector potential a^, as was the case with the tensor 
of Eqs. (44) and (46). The achievement of this gauge- 
invariancc property of the tensor [Eq. (47)] was the 
motivation [28] for taking it as the stress-energy tensor 
(on equations of motion). It remains, however, not sym- 
metric. Substituting Eq. (42) into Eq. (47) results in 



e"" = (i - £ S ) ,f A r a - ,£00 (r& e o ) # MA r ; 



(/^a)^ A r A + TT {-\ (1 - %) 



(50) 



In an empty Minkowski space = 0), the possibility 
of symmetrization of the canonical stress-energy tensor 
by adding a full derivative is provided by the conservation 
of 50(3, 1) generators. (For details, we refer the reader 
to Rcf. [26] and references therein.) In our case, there is 
a lack of isotropy due to the external field which implies 
that only a subset of the Lorentz generators are conserved 
(see the Sees. IIIC and IV A). This fact, therefore, pre- 
vents us from obtaining an equivalently symmetrized ver- 
sion of Eq. (44), and some dramatic differences arise in 
comparison with the case of an empty vacuum. The rela- 
tion obeyed by Eq. (50) to substitute for the symmctricity 
property is /0 = T /, where and / are the matrices 
and f^, while T is the transposed matrix. When 
there is no Lorentz breaking, the matrix commutes 
with the matrix /, because in this case the stress-energy 
tensor is built only of the field tensor and the unit metric 
tensor, its explicit dependence on the coordinate vector 
not being admitted. In this case, the symmetricity of 
is in agreement with the above relation. 

The stress-energy tensor 0**", as well as T^ v , satisfies 
the continuity equation with respect to the first index on 
equations of motion, 



d^ v = 0, 



(51) 



with their difference satisfying the same property 
d^dxK ^ v = 0. Its components, explicitly, are 



00 = Id ■ e + \h ■ b + a V • d 



(52) 



lJ = - d l e J - b l b? + -5 lJ (e-d- 



dod* - (V x hy 



bh) 



(55) 



where the electric and magnetic induction vectors are 
defined in Eq. (19). Explicitly, 



and 



d = (1 - £ 5 ) e + £<50 (e B) B for > 0, 

d = (l-Zs)e + £ m (e-E)E for 5 < 0, 

h = (l-£,$)b-£,$z(b-B)B for 5 > 0, 

h=(l- £ s ) b - £ 0l5 (b ■ E)E for 5 < 0. 



(56) 



(57) 



Let us consider the case in which the equations of 
motion are fulfilled, i.e., when the current j 11 = van- 
ishes identically. By integrating the continuity equation 
[Eq. (51)] with v = over a final spatial volume V and 
defining the energy in this volume as 



9% 



we get 



J d 3 xQ 00 (x) = J d?x(hi-e + \h-b 



dx° 



f Qiodai — f (e x h)id(Ti 



(58) 



(59) 



where the integral in the right-hand side is run over the 
surface S surrounding the volume V. Therefore, the 
Poynting vector zO = (e x h) % accounts for the energy 
per unit of time, per unit area, transported by the small 
electromagnetic waves. It is parallel to the group veloc- 
ity of eigenmodes and to their center-of-mass velocity, as 
described in Appendix D. In the infinite-volume limit, 
V= oo, and under the assumption that the fields e,h fall 
off at the spatial infinity, we find that the energy inside 
the infinite volume 3^° = does not depend on time: 







(60) 



By integrating the continuity equation [Eq. (51)] with 
u = j over a final spatial volume V, we find that the 
momentum in this volume. 



d 3 x e 0l (x) 



d 3 x (d x b) j , (61) 



satisfies the equation 



d& 3 v 
dxo 



Q l3 da u 



(62) 



0" 



[dxb] 1 ~a l V -d 7 



e M = [e x hy 



<<() 



d cf - (V x hy 



(53) 



(54) 



which indicates that the total momentum contained in- 
side the infinite volume ^ = ^00 conserves, 



dx° 



0. 



(63) 
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under the same assumption that the fields decrease at 
spatial infinity. Note that the momentum density ~ d x b 
and the Poynting vector ~exft describe different quan- 
tities, which does not take place in an empty space-time. 
Observe that, on the equation of motion, the structure 
of Eqs. (52)-(55) does not differ from the case of light 
propagation in an anisotropic material [33] . 2 

In the general case in which the current does 
not vanish, the spatial integral of Eqs. (52)-(53) de- 
fines the translation generators. For further convenience, 
we express the latter in terms of the momentum n l = 
dJtf /d (docii), canonically conjugated to the field a, taken 
as a canonical coordinate. It coincides with the elec- 
tric induction, defined in Eq. (19), d = —iv. We invert 
Eq. (15) so that the electric field of the wave can be ex- 
pressed as ei = s^j dj, where e" 1 must be understood as 
the inverse of the tensors given in Eq. (39) 



1 



l-£ 5 
1 

1 -£* 



1 — £5 + 23^00 
1 - £3 - 2$£, ss 



-BiBj 



EiEj 



3>0, 



5< 0. 



(64) 

With these details in mind, the translation generators 
turn out to be 

& = - J d 3 x {n x b - a (V • tt)} , 



(65) 

^° = Jd 3 x j^er.V,- + \h^b 3 - a (V • tt) J . 

Some comments are in order. First of all, we point out 
that these generators and their respective translational 
chargess Eqs. (61) and (58), differ from each other in the 
terms which are proportional to Gauss's law, V • 7r = 0. 
Such terms are intrinsically associated with the con- 
strained Hamiltonian formalism [19] (see the next subsec- 
tion). Correspondingly, wc shall show that the expression 
involved in Eq. (65) canonically realizes the space-time 
translations, at least when acting on a phase space de- 
fined by constraints associated with the gauge symmetry. 



B. Gauge fixing and Dirac brackets 

The local approximation [Eq. (43)] of our effective La- 
grangian Jzf does not depend on the velocity dao/dx°, so 
that the related momentum vanishes identically. Obvi- 
ously, this leads us to introduce 







(66) 



as a "primary constraint." Note that the symbol « must 
be understood as "weak equality"; i.e., the constraints 
cannot be assumed to equal zero until the Poisson bracket 
between two arbitrary functionals £2 and Q of the field 
variables (0^,71^) is calculated: 



d 6 x 



5£ SQ 



SB SQ 



Sa^(x) St:^{x) 5'K f j_(x) Sa^^x) 



(67) 



According to the Dirac algorithm, ipi must be imple- 
mented within the canonical Hamiltonian 



(68) 



by means of a Lagrangian multiplier so that the "to- 
tal" Hamiltonian turns out to be 



,^=3»°+ / d 3 x tf{x)ir°{x). 



(69) 



In this context, the equation of motion for J2 reads 
dj& dJ2 

and, in particular, for 21 = ttq, one has 

^ = K^}=v... 



(70) 



(71) 



The constraint ipi should hold at all time. In conse- 
quence, Gauss's law 



ip 2 = V • 7T W 0, 



(72) 



which is one of the field equations, j° = V d = 0, arises as 
a "secondary constraint." However, the latter is already 
present in J#c = in the form a V • ir [see Eq. (65)] . 
Therefore, ao can be considered as a Lagrange multiplier 
and thus an arbitrary function of x. Its equation of mo- 
tion implies that 



dao 
d^ 



{a Q ,J^} = tf(x). 



(73) 



The string of constraints stops here because Gauss's 
law commutes with the Hamiltonian. Moreover, our pri- 
mary and secondary constraint are " first class" with a 
vanishing Poisson bracket: 



{7r (x),V-7r(x')} = 0. 



(74) 



The remaining algebra of the constraints and the Hamil- 
tonian is given by 



{Tro(x),Tr (x')} = 0, {V • ir(x), V • tt(x')} = 0, 

{7T ,^c} = -V-7T, {V-7T,jr C } = 0. 



(75) 



2 Remember that we restrict ourselves to the low- frequency, low- 
momentum limit in the present subsection and in the rest of the 
article. 



Observe that the Lagrangian multipliers & and ao trans- 
fer an arbitrariness to the Hamiltonian [Eq. (69)]. As 
a consequence, we are forced to deal with a phase space 
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plagued by nonphysical degrees of freedom. This problem 
is closely associated with the gauge invariance property 
and is formally removed by imposing two gauge-fixing 
conditions. Because of this fact, the two existing multi- 
pliers are rendered to precise dependences on the physical 
fields, and can eventually be removed from the theory. 

A suitable set of gauge conditions can be found by 
solving Gauss's law [Eq. (72)] with respect to a : 



1 



a 



VeV 



do (Vieyaj) 



(76) 



where VeV = V;£yVj. Guided by this result, we are led 
to choose a generalized version of the Coulomb gauge as 
the third constraint of the theory 



<P3 = ^iEijdj 



0. 



(77) 



The consistency consequence of this gauge condition 
can be found by Poisson-commuting the Hamiltonian 
[Eq. (69)] with ip^. However, it can be read off directly 
from Eq. (76) and promotes the last constraint, 



ifii — a pa 0. 



(78) 



Since this is found to be stationary as well, Eq. (73) pro- 
vides a vanishing value of the Lagrangian multiplier ff, 
and there are no further constraints. 

The accessibility of these gauge conditions can be 
checked by using the gauge function 



A 



(79) 



In fact, for any value of a, a , the gauge-transformed 
fields 

afe) = a t {x) + ViA, a' (x) = a (x) + d A (80) 

obey the gauge conditions of Eqs. (77) and (78). Explic- 
itly, 3 



V.^ija'j = Viffyflj + V l e ij V j A = 0, 
V • 7r' = -ViSijVja'o = ViEijVjaQ + V, = 0. 



Certainly the set {<fi} is second class with a character- 
istic matrix C a p{x,x') = {(p a (x), ipp(x')} given by 



Cat 



1 

\7 x e\7 x 

-V x eV x 
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6 {3 \x-x'). (83) 



Since C a p(x,x') is regular by construction, we can also 
write down its inverse 



/ 



a 



a0 











-1\ 



1 








S (3) (x-x'). (84) 



J 



With the help of the latter we introduce the Dirac 
bracket: 

{3{x), C(*')}. = {-W. C(*')} - / d 3 y {£(x),ip a (y)} 

|rf 3 zC- /3 1 ( J /,z){ W (z),Q(a ; ')} (85) 

In this context, the fundamental bracket of the theory 
can be calculated straightforwardly and reads 



{ai(x),irj(x')}^ = Ujix.x') 



where 

tij {x , X ) = 



x x 



(86) 



(87) 



is a projector- valued distribution which fulfills the rela- 
tion 

J d 3 ytu(x,y)ti j (y,x') = tij(x,x'). 



(81) 



Note, in addition, that the constraints {ifi} defined above 
restrict the original phase space of the theory to a four- 
dimensional hypersurface, 



Q = {(a M ,7r M ) | tpi w 0, i= 1...4} 



(82) 



This, however, is not symmetric tij(x,x') ^ tji(x,x') A 
and, depending on the field to be projected, one must 
contract one index or the other. For instance, let us de- 
compose the canonical field into two mutually orthogonal 
pieces: a(x,t) = a T (x,t) +a e (x,t). Here the transversal 
component is obtained by contracting the field a(x,t) 
with the index of ty ;(x,x'), which is provided by the op- 
tical tensor ey. Thus, 



in which the time evolution of two physical degrees of 
freedom takes place. 



aj{x,t) 



J d s x'tij(x,x')aj(x',t) 



'VeV 



VfcEfcj 



aj(x,t) (89) 



3 It is worth observing at this point that the electric field associ- 
ated with the small-amplitude waves e = Vao — doa is a gauge- 
invariant quantity. As a consequence, the canonical momenta 
Hi = —Eijej are invariant as well. So, under the gauge transfor- 
mation of Eq. (80), 7r' = it. 



4 For a vanishing external field, this reduces to the symmet- 
ric transversal projector associated with the Coulomb gauge 
tf™\x,x') = [8 i:j - ViVj/V 2 } S^ix-x'). 
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Accordingly, the longitudinal components of a turn out 
to be a\{x,t) = J d 3 x' Cij(x,x')a,j{x' ,t), where 

^■(x,x')^Vf^V^ fcj( 5( 3 )(x-x') (90) 

is the nonsymmetric longitudinal projector. Likewise, we 
decompose 7r = ir T However, to extract its transver- 
sal and longitudinal elements, the canonical momentum 
must be contracted with the index provided by the gra- 
dient operator: 

nj(x,t) = J d 3 x'tij(x,x')iTi(x',t) 



1 



1 VeV 



7Ti(x,t). (91) 



Due to the gauge-fixing condition [Eq. (77)] and 
Gauss's law [Eq. (72)], the longitudinal compo- 
nents a\ = J d 3 x'Cij(x,x')aj(x',t) and i^{x,t) = 
/ d 3 x' Cij(x,x')iTi(x' , t) vanish identically. Keeping these 
details in mind, the fundamental Dirac bracket of our 
system [Eq. (86)] acquires the following structure: 



{a[(x),^J(x')}, = *«(*,*') 



(92) 



On the other hand, the Dirac bracket of clq or ttq with an 
arbitrary functional J2 vanishes identically by construc- 
tion: 

{4(x), ao(»')}, = 0, {£(z),M*)}. = 0- (93) 

Because of this, the sector (do, ttq) can be formally elim- 
inated from the phase space, and the theory is fully de- 
scribed in terms of (a T ,7r T ). Moreover, with being a 
generic functional one has 



■-P2 



<P3 



{^(x),V-tt(x')}* = {^(x),V l e lJ a 3 (x')}^0. (94) 

Thus, both the second and third constraints [Eqs. (72) 
and (77)] are no longer "weak" equalities, and instead 
they can be used as "strong" equations. The latter termi- 
nology is conceptually equivalent to replacing the symbol 
w by = in the set of constraints {fi}- So, once the Dirac 
brackets [Eq. (85)] are constructed, they can be set to 
zero everywhere. 

Making use of Eq. (85) and Eqs. (88)-(91), we find 

{^°,a;(z)h = d°a t (x), {&°, n (x)}, = d^x), 

{&>, ai (x)}^ = V ai (x), {&>, 7ri(x)} t = Vm(x) 

which are the well-known time and spatial transforma- 
tion properties of fields. It follows that for any poly- 
nomial functional £2 of a and ir that does not depend 
explicitly on x, one has 



{&>i*,£(x)}« = d"£(x) 



(95) 



To conclude this subsection, we derive the modified 
Maxwell equations. Whatever the nature of the electro- 
magnetic background (5" > or 3" < 0), the Hamiltonian 
equation of motion for tt becomes Ampere's law: 



f/7T 



V x 



6b 



-V xh. 



(96) 



Together with the constraint if 2 = (Gauss's law), they 
make the second pair of Maxwell equations. The Hamil- 
tonian equation for 6, 



db 



for 



= -V x e, 



(97) 



where e = — dyd, becomes Faraday's equation. It is ful- 
filled as an identity. Together with another identity, 
Gauss's law for magnetism, V • b = 0, which is not a 
Hamilton equation of motion, they make the first pair of 
Maxwell equations. 

We stress that the procedure developed in this subsec- 
tion is also applicable to any other linear approximation 
of electrodynamics in which the optical tensors depend 
on neither the space-time coordinates nor the derivates 
with respect to the latter. Observe that it is even suitable 
to describe the situation in which there exists a certain 
biaxiality associated with the crossed fields configuration; 
i.e., where the external field invariants vanish identically, 
3 = 6 = 0. 



C. Generators of rotations and Lorentz 
transformations 

In this subsection we first obtain the conserved gen- 
erators associated with the Amputated Lorentz Group 
50^(3, 1). The Noether theorem for infinitesimal trans- 
formations from the SOa(3, 1) group over the field a^(x) 
that leave the action F in Eq. (2) invariant reads 



0, 



(98) 



where Jz? is the quadratic Lagrangian of Eq. (17), T fi v 
is the canonical stress-energy tensor of Eq. (44), and the 
transformation laws are 



(99) 



with tu al3 determined by Eq. (27) and with the vectorial 
representation of the Lie algebra generator of 5*0(3, 1) 



P) ■ <T =i[8 P pVc ur -S p a r,p«). (100) 



With these details in mind, the Noether conservation 
equation [Eq. (98)] may be written as 



()„ 



= 0. 



(101) 



13 



Correspondingly, the conserved currents associated with 
the SOa{3, 1) symmetry are given by 



J^l&apJ^ and > = ^#^J^. 



Their respective continuity equations read 
1 

2' 

1 ^ 



-& a pd^ = 0, 



0. 



(102) 

(103) 
(104) 



To provide the fulfillment of Eqs. (101) and (103), and 
(104), it is sufficient to define the in generally noncon- 
served current J^ a @ in the form 



j/xa/3 _ _ 



T^x a 



<9(<9 M a„) 



(105) 



which imitates the elements associated with 50(3, 1) in- 
variance. For further convenience, we express J^ a P in 
terms of Q^ v [Eq. (50)]. To this end, we substitute 
Eq. (48) and make use of the identity 



as? 



dj£_ 
<9(3 M a A ) 



Consequently, 

jm"/3 = -iQ^ x " (y^)^ = x a Q^ - 3^0*". (106) 

Observe that the nonconservation of J^ a @ is intrinsically 
related to the nonsymmetric feature of the energy mo- 
mentum tensor [Eq. (50)]. In fact, 



(107) 



Since the continuity equations [Eqs. (103) and (104)] in- 
volve the projection of Eq. (107) onto the external field 
tensors, we end up with the following identities: 



& a pQ afi = and ^ a B& af3 = 0. 



(108) 



Now, the respective spatial integrals of the time com- 
ponents of the currents contained in Eq. (102) provide 
the conserved charges 



1 



with 



dx° 



and Sf = -JV £ 



1 ~ 

2' 



and 



93 



= 0. 



(109) 



(110) 



These scalars involve a second-rank tensor whose struc- 
ture extends the known representation of the Lorentz 
generator 



d 3 xJ aflu (x,t) 

d 3 x (x^Q 0l/ - x»e ^) (111) 



to the case of the violated Lorentz invariance under con- 
sideration. Correspondingly, we can define both the pho- 
ton angular momentum ^ = (J f23 , ^X 31 , ^ 12 ) and the 
photon boost generator JT = \j w , c / 20 , f 3 ®). 

Considering the prescription above we can express the 
conserved charges of SOa{3, 1) in the most general case 
of a constant and homogeneous external field as 



& = b- s -E-jtr, 
<i = BJtr+E a f. 



(112) 



Thus, as soon as the effects of the vacuum polarization 
tensor are considered the number of Lorentz generators 
is reduced from 6 to 2. Therefore the vacuum symmetry 
group ISOa(3, 1) has dimension 2. 

The explicit structure of ^ and Jfc in terms of the 
canonical fields (a,w) follows from Eqs. (Ill) and (47), 
and reads 

J = - J d 3 x{x x[n xb}- {x x o)(V ■ tt)} , (113) 

Jf = -x°& + J d 3 x jx Q- , + \h^b 

- aoV-7r)}, (114) 

where b = V x a. Note how the secondary constraint of 
our problem is implemented in both generators through 
terms proportional to ~ V • it. In the phase subspace 
defined in Eq. (82), such terms vanish identically, and 
the resulting expression of ^ coincides with the angular 
momentum of light in an optical medium [33] . 

Now, using the definition of the Dirac bracket 
[Eq. (85)] and equipped with the "Lorentz-like" gener- 
ators and J^i in Eqs. (113) and (114), we are able to 
express 

{^^(x)}* = (x x V)iflj(x) + €ijiai(x) 

-e»fcm VgV o,i{x) - e ikl -y-y — a ; (x(j;i5) 
{Ji^jix)}* = {xx V) l TT j (x) + eiji7Ti(x), (116) 
{J£j, a 3 (x)}* = x l d°aj(x) - x°\7 l a J (x) - ^ 7 ^ kl ^ 1 

x [ Xl d a k {x)] , (117) 
{J^,7Tj-(x)}* = -x°Vi7rj - e l]k ^ n (V x a) m 

+Xie jkl Hi^k (V x a) m (118) 

These expressions deserve some comments. The first 
pair of brackets realizes the infinitesimal rotation on the 
canonical variables a and 7r, respectively. Note that, in 
contrast to the momentum 7r, the field a transforms as a 
vector up to terms associated with the the gauge-fixing 
condition [Eq. (77)]. A similar statement applies to the 
boost transformation properties contained in Eqs. (117) 
and (118). We find it convenient to emphasize that all 
these brackets guarantee that the infinitesimal canoni- 
cal transformations induced by the generating functions 
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, J£i do not lead out from the constrained phase sub- 
space ofEq. (82). Indeed, 



^m^mjl^j (»e) }* 



0, V,{ 

: 0, V J {X i ,'n j {x)} i 



0, (119) 
0. (120) 



While the first column realizes ip^ = [Eq. (77)], the sec- 
ond one verifies (f2 ~ [Eq. (72)]. It is worth observing 
at this point that in the limit in which the external field 



vanishes, £; 



Sij , the gauge condition is reduced 



to the standard Coulomb gauge (V • a = 0). Due to this 
fact, the rotation transformation property of the gauge 
field [Eq. (115)] becomes similar to the one associated 
with the canonical momentum [Eq. (116)]. 

The transformation properties of the gauge field a(x) 
in Eqs. (115) and (117) are not very helpful as they 
stand. For computing more cumbersome brackets in- 
volving gauge-invariant quantities, it is more convenient 
to represent the a-containing part of Eqs. (115)-(118) in 
terms of the magnetic field h and magnetic induction b, 
respectively. To this end we apply Vx to Eqs. (115) and 
(117). Then 

{fu bj(x)}* = (xx V)ibj(x) + e ijk b k (x), (121) 
{Ji^jix)}^ = (x x V)i7Tj(a;) + e ijk ir k (x), (122) 



{jei,bj(x)}* 
{Jfi, ^(x)}* 



-e ijk e k + x*(V x e)j - x V 4 6/,123) 
-e ijk h k + x 4 (V x h)j - .t°V 4 t(^24) 



These brackets constitute the starting point for deter- 
mining the effects induced by the vacuum polarization 
within Lorentz algebra. Note the remarkable symmetry 
under the interchange 7r ^ b of the first pair of these 
equations. The second pair, however, is invariant under 
the simultaneous replacement 7r ^ b and e ^ h. 



IV. DIRAC COMMUTATORS AND 
EQUATIONS OF MOTION FOR GENERATORS 

The goal of this section is to determine the equations 
of motion of the generators of the Lorentz rotations and 
to establish the Dirac commutation relations distorting 
the Lie algebra of the Poincarc group. By contracting 
these relations with the field tensor and with its dual, 
the closed algebra of the vacuum symmetry subgroup is 
obtained. As with other algebraic relations, the vacuum 
invariance holds in the physical subspace of the phase 
space specialized by constraints, where the evolution of 
the physical degrees of freedom takes place. 



A. Equations of motion for the angular momentum 
and for the boost generator 

The vacuum polarization tensor provides an effective 
coupling between photons and the external field # '. In 
order to explore the role of this quantity within LSB, we 



first consider the equation of motion associated with the 
total angular momentum of the electromagnetic waves: 



dx° dx° 



1 ) H 



(125) 



The first term on the right-hand side vanishes identically, 
since ^ does not depend explicitly on time (x° = t), 
whereas the Dirac bracket provides the total torque ex- 
erted over the photon field. In our context (see details in 
Appendix A 1), this can be expressed as 



{</<^°}* = / d 3 x[ir xe + hxb] 



(126) 



In the rotation-invariant case-say, a vacuum or isotropic 
material-whcn the dielectric permeability is a unit ten- 
sor, Eq. (126) disappears in correspondence with the mo- 
mentum conservation. In this case, 7r||e, h\\b, and their 
vector products are zero. As a consequence, the equation 
of motion for the photon angular momentum can be writ- 
ten in terms of the spatial part of the energy momentum 
tensor [Eq. (55)]: 



dx° 



1 

r 



ijk / r/ViB' 4 (-)'' ■' :. 



(127) 



On the contrary, the equation of motion associated 
with the photon boost generator is given by 



dJXT 

dx° 



dJtr 



(128) 



Because of the explicit dependence of on time [see 
Eq. (114)], the first term on the right-hand side con- 
tributes to the equation of motion dJf /dx° = — & . The 
Dirac bracket involved in this expression is calculated in 
Appendix A 2 and reads 



{jfT,^}^ = J d 3 x exh. 



(129) 



Combining these details, Eq. (128) acquires the following 
structure 



dJT_ 
dx° 



J d 3 x[(e x h)* + (tt x bf] 

J d 3 x[Q l ° - 6° l ] (130) 



where 0*° and 0° l are the Poynting vector and the den- 
sity of momentum, respectively. 

Both Eq. (127) and Eq. (129) can be embedded in a 
compact four-dimensional expression 



dx° 



(131) 



in which the antisymmetric tensor 1^ is connected to 
the stress-energy tensor by means of the following expres- 
sion 



d 3 3 



'(*), 



(132) 
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The structure of Eq. (131) with Eq. (132) is somewhat 
expected: the spatial integration of Eq. (107) reproduces 
the same equation for the Lorentz-like generators up to a 
surface integral ~ <f s doi J lQ/3 , which vanishes identically 
when the rapid falloff of the canonical fields at spatial 
infinity (S — > oo) is provided. 

Observe that the projections of Eq. (131) onto the ex- 
ternal field tensors are consistent with the conservation 
law of Sf and S? [Eq. (110)]. In fact, due to Eq. (108), we 
obtain that 

Ij^T^O and i#^T^ = 0. (133) 

We shall see that these identities, together with the equa- 
tion of motion of ^^ u , turn out to be very convenient 
to evaluate the vacuum polarization effects on "Poincarc- 
like" algebra. 



B. Distorted Poincare algebraic relations 

Let us consider Eq. (95) with replaced by 8 0iy . After 
an integration over x we obtain 



\P"\= f d 3 xd p O Ql/ {x,t). 



(134) 



Using the divergence theorem and assuming that the 
canonical fields a and tt vanish sufficiently rapidly at in- 
finity, wc find 

J d 3 xd l e 0,J (x,t) =0. (135) 

Therefore, 

{& ,i ,& ,j } t: = and =°- ( 136 ) 

Because of the antisymmetry of the Dirac bracket, we 
have that ^ g0 } t vanishes identically as well. Hav- 

ing these aspects in mind, we can write Eq. (134) as 

{&n* t &"'} m =0. (137) 
We also apply Eq. (95) to the following case: 

J 0Xa } r = x x {^, e ^ - x G e 0A }^ 

= x x d p O a,J - x CT ^e 0A (138) 

where the integrand in Eq. (Ill) has been considered. 
Thanks to the identity x a d (S Q° T = d^(x a Q Or )-r] 0a e OT , 
we are able to express Eq. (138) as 

{^,j ox<t }^ = v pa e ax -T/^e ' 7 - ^ J 0A<7 .(139) 

Integrating the expressions above over x, we find 
| ^ Act } = r^® gfi x — r\ pX 

d 3 xd p f )XlJ . (140) 



Because of the rapid vanishing of the field at infinity, we 
obtain that, for \i = i, J d 3 xd t J 0X(7 vanishes identically. 
However, if /j = 0, the last integral can be written as 
d° J d 3 xJ 0Xa =d Q f Xa . With these details in mind, and 
using Eq. (131), we end up with 

{^W / x °}^ = jf°@> x - jf x &° - if Q 1 X(J , (141) 

where 1 Xa is specified in Eq. (132). The bracket above 
reproduces the commutators associated with ISO(3, 1) 
Lie algebra up to a term manifesting LSB. 

Our analysis on the Poincarc-like algebra is to be com- 
pleted by deriving the Lorentz-like algebra, which in- 
eludes J^}^, and {jr%Jf^. A 

detailed derivation of these Dirac brackets can be found 
in Appendix B. In particular, we have obtained that 



{J 1 ,^ }^ = e ijk JtT k + -e llm I d 3 x xV m , ^ 42 ^ 



with t ij being the spatial part of the tensorial density 
T pv defined in Eq. (132). Naturally, the brackets above 
can be combined in a four-dimensional expression, 

l^/A^ J> pa } = r\ pp J> va —jf p + i]' Jfi J? pl> 

-rf V JP* - J d 3 x(x"^° - X^Tf Q )T pa 

d 3 x{x p if* - x a r] p0 )T^ (143) 

where the the right-hand side reproduces the standard 
well-known result of the the Lorcntz algebra up to 
terms that contain the antisymmetric part of the energy- 
momentum tensor. Note that in spite of its nonconser- 
vation, the spatial rotation generators retain the stan- 
dard angular momentum algebra according to the first 
line in Eq. (142). This is a usual situation in classical 
and quantum mechanics with broken symmetries under 
transformations that do not touch the time variable. It 
cannot be the same with the whole of the Poincare alge- 
bra, because the Hamiltonian itself is one of its genera- 
tors, whereas its commutators with other members of the 
algebra that do not conserve cannot help being affected. 
That is why the violation of the algebra in the second and 
third lines of Eq. (142) is not unexpected. This result in- 
dicates, in addition, that the Lorentz invariance in pure 
Coulomb-gauge V • a — Maxwell theory is established 
only when the Lorentz generators are conserved. We con- 
clude our analysis by pointing out that the brackets in 
Eqs. (141) and (143) coincide with algebra obtained in 
Ref. [21] within the context of the free field quantization 
in a noncovariant gauge. 



16 



C. Algebra of the symmetry subgroup 



In order to pursue our research, we proceed to contract 
Eq. (141) with \^\ a and \&\ a - As a consequence, 

{Sf,^h = ^^-J^ij^T^, (144) 

jV,^ 1 } = /""^-j]" ^^, (145) 
I J * 2 

where the definitions of and Sf have been used. Accord- 
ing to Eq. (133), the last terms on the righ-hand sides of 
these brackets vanish identically. Therefore, 

.^.'Pv, = (146) 

The remaining Dirac bracket involving and §f can be 
determined by projecting Eq. (143) twice onto the exter- 
nal field tensors. This procedure generates the following 
brackets: 



u P a 



(147) 



Both J^a^ A and A = 2^ + i^jr, a 

are symmetric tensors. Hence, their contractions with 
the antisymmetric tensor vanish identically, and 

{^,5?}* = = 0. The latter result is expected 

because it comes out from Eq. (67) that {J2, £2} = 0, 
with J2 being a generic function of the canonical vari- 
ables. Furthermore, by considering Eq. (133), we can 
claim that the right-hand side of the last Dirac bracket 
vanishes identically as well. 

We can then summarize the Lie algebra of ISOa(3, 1) 
as follows: 



= o, 



= 0, 



(148) 



with the external field tensors & and & playing the roles 
of group structure constants. Certainly the translation 
generators induce an Abelian-invariant subalgebra which 
defines the nonsemisimple structure of 50^(3, 1). More- 
over, the Casimir invariants of our problem are given as 



0> 2 = Z 1 +Z 2 , 

Zi = and Zi = - 

23 



23" ' 



(149) 



It is worth mentioning at this point that the scalars in- 
volved in Eq. (9) are, therefore, maps of the invariants 



above: & 2 h-» k 2 , Zi H> zi, Z 2 H- z 2 . Moreover, in the 
special frame where the field is purely magnetic ($ > 0) 
or purely electric ($ < 0), Eq. (143) expands into 



0, 



where Eq. (112) has been used. Each column in this 
set of commutators manifests a subalgebra: the first 
one corresponds to the two-dimensional Euclidean group 
ISO(2), whereas the second one corresponds to (1 + 
1)— dimensional pseudo-Euclidean group ISO(l, 1). The 
latter groups arc associated with the transverse and pscu- 
doparallel planes with respect to the B(E) direction. 
Therefore, the symmetry subgroup SOa{?>, 1), down to 
which the Poincare group is broken due to the presence 
of an external field, reduces in the reference frame-where 
that field is purely magnetic or clcctric-to the direct 
product of ISO(2) and ISO(l, 1). Besides, we want to 
remark that as long as a photon propagates transverse 
to the external field, the square of the Pauli-Lubanski 
operator zu^ = 1 / '2e MAo ' e ^ \ a & e is no longer a Casimir 
invariant. This fact reflects the underlying difference be- 
tween the vacuum in an external field & and the case 
of an empty space-time in which the SO (3, 1) symmetry 
is preserved and all particles are classified according to 
the spin and hclicity representations encoded in the w 2 
eigenvalues. 



V. MAGNETIC MOMENT OF SMALL 
ELECTROMAGNETIC PERTURBATIONS OF 
THE VACUUM 

In this section we analyze some consequences associ- 
ated with the equation of motion for the photon angular 
momentum [Eq. (125) and (126)]: 



dx< 



\=1 with 1 = J d 3 x(n x e + h x b). (150) 



In the special case where the external field is magneticlike 
(3 > 0), the explicit substitution of n and h [Eqs. (56) 
and (57)] into T [Eq. (150)] allows us to express the latter 
in a rather meaningful form: 



CM = - / d 3 x 

2 ./ I e i e 



1 = 2M x B 

£0® 



(ir • B) 7T + % 5 (b-B)b 



(151) 



where Eq. (40) has been used. The expression above 
mimics the torque exerted by the external field on the 
magnetic dipole iW. The torque 1 [Eq. (151)] vanishes 
when projected onto B, so Eq. (150) implies that the par- 
allel component of the photon angular momentum, &\ , 
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is a constant of motion. On the contrary, the projec- 
tion (helicity) f) ~ ^ ■ & of the angular momentum of a 
photon onto its canonical momentum [Eq. (65)] 5 is not a 
conserved quantity unless & turns out to be parallel to 
the external field. 

The magnetic moment fW is a feature of the small elec- 
tromagnetic perturbations of the vacuum (photons in the 
first place), when the former interact with it through 
virtual electron-positron pairs. This interaction makes 
the photon behave like a magnetic dipole. Note that M 
[Eq. (151)] is a gauge-invariant quantity, also orthogonal 
to the photon canonical momentum, & ■ M = 0. More- 
over, it has two components in correspondence with the 
cylindrical symmetry imposed by the external magnetic 
field. One of them, 



M± = ^ f d 3 x j£ 00 ^-^7r ± + £ 55 {b -B)b ± 



is perpendicular to B; whereas the remaining one, 

flf H = \ J d 3 x | ilee !^tt|| + £ SS (b ■ B) b I (153) 

is parallel to the axis in which the external field lies, and 
therefore invariant under rotation about the magnetic 
field direction. This, however, does not contribute to 
T [Eq. (151)]. Hence, it does not play any role within 
the equation of motion of the photon angular momentum 
[Eq. (150)]. Nevertheless, contributes to the effective 
Hamiltonian [Eq. (65), H = and thus to the photon 
energy. Note that with the use of Eq. (40), the latter can 
be conveniently written as 



H = 



1 



T/'J 



-U2 



MB. 



(154) 



In contrast to the perpendicular component of M 
neither remains invariant under a rotation around the 
B direction nor contributes to the photon energy since 
it is projected out from the scalar product involved in 
Eq. (154). However, it turns out to be a clear manifes- 
tation of LSB, since it specifies, by means of Eq. (150), 
that not all components of the photon angular momen- 
tum are conserved quantities. Furthermore, it follows 
from Eqs. (40) and (153) that fMj| > behaves param- 
agnctically. By contrast, it is not possible to establish 
a definite magnetic behavior in fM± because it contains 
terms which mix not only 7T|| and 7rj_, but also b± and 6|| . 

Let us consider the case in which the external mag- 
netic field is asymptotically large, b = |B|/B C — > oo. In 
this limit, the basic entities contained in Eq. (154) [see 



5 Recall that this is the direction of the momentum flux and the 
wave vector (see Appendix C), not of the energy flux, whose 
direction coincides with the Poynting vector and with that of 
the group velocity. 



Eq. (41)] are written in Rcf. [28], referring to £ as the 
one-loop term of the Eulcr-Hcisenberg Lagrangian [34- 
37]: 



2S£ C 



£ s «-m & , 



OL Oi 



(155) 



Observe that for a magnetic field with 10 < b <C 37r/a 
one can treat e±_ = fij 1 ~ 1 and therefore tt ~ — e — 
£gg (e-B)B. The resulting effective Hamiltonian H reads 



H 



d S X [ -7T 

2 



1 -b 2 
2 



MB. 



(156) 



In this approximation, the second term in Eq. (151) de- 
creases as 1/b, and thus contributes in if as a small con- 
stant. Hence, it can be disregarded in comparison with 
the term provided by the first term of Eq. (151), which 
turns out to be a linear function on the external field 
strength. As a consequence, the magnetic dipole acquires 
the following structure: 



M^g—S, 

2m 



(157) 



with £ = a/3ir being a sort of Lande factor, whereas 









L 


i / 


d 3 







(158) 



It is worth mentioning that «S is only determined by the 
electric induction vector associated with the small elec- 
tromagnetic waves. Thereby M can be interpreted as a 
magnetic moment, with S playing the role of the "spin" of 
the small electromagnetic waves. This terminology, how- 
ever, is used just to establish an analogy with the case of 
the electron magnetic moment. In contrast to any mas- 
sive particle, a photon lacks a rest frame, and thus one 
cannot define a spin for it. Besides, S does not fulfill the 
standard Dirac bracket of the angular momentum, i.e., 
the first bracket in Eq. (142). 

We continue our research by considering the electric 
field of each cigenmodc as a monochromatic plane wave, 
so that 



\eW(k) 



COS[LO\X 



(I 



k ■ x] 



(159) 



Here and ui\(k) are the amplitude and frequency 

of mode A, respectively. Note that the shape of (k) 
can be found below Eq. (13). Observe, in addition, that 
only mode 2 has an electric field parallel to B. As a con- 
sequence, S becomes physically relevant for the second 
polarization mode. Inserting the expression above into 
Eq. (158), we find 



d A x 



k±'E 



(2)2 



m|fc| 



-5cos 2 [a;2X 



(160) 
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where the leading term has been withheld so that a/ 2 ' sa 
|fc and s = /\e^\ w n x (ny xnj. Equipped with 
these approximations, the time average of fW reads 



(161) 



Here, V denotes the volume over which the integral con- 
tained in S is considered. It is worth observing at this 
point that fW is proportional to the average energy asso- 
ciated with the second propagating mode: ~ |£q . 

We should also mention at this point that the notion of 
the photon magnetic moment was introduced in Ref. [23] 
and discussed in Refs. [24] and [25] . In those works, it was 
defined as contributing to the photon energy as a function 
of both its momentum and the external field strength. 
Moreover, it was conceptually analyzed in different en- 
ergy regimes and magnetic field contexts of the photon 
dispersion curve. Only in Ref. [25] were the asymptotic 
conditions investigated in the present work considered. 
The connection between the respective photon magnetic 
moment and iW in Eq. (161) can only be established un- 
der the second quantization of our problem, which re- 
quires us to substitute the Dirac brackets with standard 
commutator relations. 6 In such a case, the frequencies 
u>2(k) and ujs(k) are related to the respective photon en- 
ergies, and in correspondence, the mode 2 photon density 
turns out to be 3\£ = ^£q 2 ' 2 /o; 2 . Dividing (fW) by the to- 
tal number of mode 2 photons contained in the volume 
V, wc obtain 



(AO e rn x 



(162) 



where m is the electron mass and f{k±_) = fcj_/m is a 
dimensionless form factor which guarantees the gauge in- 
variance of the theory; i.e., m does not provide a photon 
rest mass m(k± — > 0) — » 0. Equation (162) coincides with 
the photon anomalous magnetic moment previously ob- 
tained by one of the authors in Ref. [25]. However, in 
that work, tn was defined as the coefficient of B when 
the dispersion curve u)2{k) is linearly approximated in 
terms of the external magnetic field. 



VI. CORRESPONDENCES WITH THE 
GENERAL LORENTZ -VIOLATING 
ELECTRODYNAMICS 

The present work is not concerned with all thinkable 
Lorentz violations associated with extensions beyond the 



6 The second quantization of the small electromagnetic field in an 
external field becomes a necessary issue as far as one wishes to 
go beyond the purely electromagnetic sector by including not 
only virtual, as here, but also free charged particles. This task, 
however, is not the issue the present work, although some needed 
building blocks arc prepared in it. 



Standard Model like Refs. [9-11]. Just the opposite; it 
is developed entirely within Standard Model and deals 
specifically with such Lorentz symmetry violations as are 
stimulated by a background electromagnetic field, mostly 
by a time- and space-independent magneticlike field (<3 
= 0, 5 > 0) within quantum electrodynamics (QED). 

Nevertheless, it makes sense to try to reduce the ap- 
proaches of Refs. [9-11] and of related studies [henceforth 
referred to as General Lorentz- Violating Electrodynam- 
ics (GLVE)] to a common denominator with the results 
of many works dealing in a rclativistic, covariant way 
with external fields, nontrivial metrics, and/or a medium, 
as was proposed by a referee of PRD. When treated in 
the framework of conventional physics, these are acting 
as Lorentz and £0(3) invariance-violating agents, and 
therefore may supply special examples to serve as models 
for verifying general constructions in GLVE. A full anal- 
ysis in this field would require a quite separate study. So 
as not to deviate too far from our principal theme, we 
are now only listing-for an external magneticlike ficld- 
the CPT-even (kf) k>j ' Xi ' coefficients, whose combinations 
are subject to measurements in various experiments in- 
tended for detecting Lorentz violations, as they follow 
from the general covariant decomposition of the polar- 
ization tensor in a nonlinear electrodynamics. We shall 
see that, contrary to the postulate accepted in GLVE, in 
our context this tensor is not double traceless, but quite 
the opposite; its trace is physically meaningful and asso- 
ciated with the diffractive properties of the anisotropic 
medium formed by the background field. Even if admit- 
ted, the case where the double trace would vanish iden- 
tically could not introduce any modification to the dis- 
torted Poincare algebra, i.e., Eqs. (141) and (143). The 
double tracelessncss condition modifies the optical ten- 
sors of the theory, but the Poincare-like generators keep 
their structure as long as they are expressed in terms of 
the canonical variables. Besides, the energy-momentum 
tensor remains nonsymmetric, also as in Ref. [9], a fact 
needed to save the nontrivial equation of motion of the 
angular momentum of light. 

We also comment on how the sensitivity achieved in 
experiments aimed at detecting possible Lorentz viola- 
tions in the vacuum might be confronted with measur- 
ing equivalent magnetic fields. We made sure that the 
accuracy available would not be sufficient to detect the 
Lorentz violation produced even by the hitherto strongest 
laboratory magnetic field in the vacuum. To detect an ef- 
fect of Lorentz violation presumably inherent in the vac- 
uum, which might be equivalent to the one stemming 
from QED with an external magnetic field on the order 
of the cosmic background (10~ 6 G), the experimentalist 
would need to achieve a sensitivity of 10~ 44 in experi- 
ments that might, besides, exclude the influence of any 
dielectric material involved in the experimental device. 
This surpasses the boldest prospects of sensitivity un- 
der present-time considerations by at least 30 orders of 
magnitude. 
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A. The components of the CPT-even tensor 

The quadratic part of the effective Lagrangian £ 
[Eqs. (17), and (18) and (43)] in an external magnetic 
or electric field can be expressed as 



(163) 



In accordance with Eq. (2), this is the same as 

Thanks to the antisymmetry of the field tensor f^, the 
coefficient tensor (kf) kX ^ v is not defined by Eq. (163) 
in a unique way: its parts, symmetric under the permu- 
tations within the first and second pair of indices, are 
left undetermined. For this reason, this tensor is to be 
understood as antisymmetric under these permutations: 



(k f ) 



XhtflV 



(164) 



Analogously, only the part of the tensor (kj?) kA '" / , which 
is symmetric under the permutation of the first and sec- 
ond pairs of indices, contributes to Eq. (163). For this 
reason, this tensor is also to be understood as symmetric 
under these permutations: 

(k f T X ^ = (k f T" kX . (165) 

By using the definition f^ v = i(k fl a l/ — k t ,a )l ) in Eq. (163), 
we see that the polarization tensor that has direct phys- 
ical meaning is connected with combinations of tensor 



(k f ) 



components in momentum space as 
1 



x [(k f T vX ^ - (k f ) k ^ x - (kfY kXv + (k f Y^ x \ 
= 2k K k x (K F ) x ^. (166) 

The properties of Eqs. (164) and (165) provide that this 
polarization tensor will be symmetric, IL^, = n„ M , as it 
should be in the vacuum with a background field [27]; 
and transverse, H^k^ = 0, as is prescribed by the gauge 
invariance. Using the relation 

= 23'(?7 a!CT <^ - V a0 S a e ) + Tf°&\&\ 
-S^&^&f + 5^ aX ^ x a - rf a & aX ,9xt 

in Eq. (43), the coefficients (kf) kX ^ v for the infrared limit 
in QED with a constant magnetic field may be chosen as 



(k f ) 



1 



2^ 

J(£ es -£^)J^ A J^. (167) 



The (anti)symmetry properties of Eqs. (164) and (165) 
are obeyed by Eq. (167). The Lorentz violation induced 
by the magnetic field is characterized by 6 — 1 = 5 compo- 
nents of the antisymmetric external field tensor sub- 
jected to one condition = 0, and by the three scalars 
£gj £ei5j £53 determined by the dynamics of the inter- 
action. 

Beyond the infrared limit, Eqs. (17) and (18) imply 
that the extension of the tensor (kf) kiiXv to include the 
infinite series of the space-time derivatives (kf) kiiXv — 
En( K r) ¥Ami '"°°^-^ reduces to multiplications-in 
the momentum space-of its separate parts by the scalar 
functions f^-, *^g 2 * 2 , and ^j^f built of the polarization 
operator eigenvalues, or to equivalent action of the corre- 
sponding (nonlocal) integral operators in the coordinate 
space. This factorization feature is not a consequence of 
Eqs. (14) and (20), since these equations hold true al- 
ready in a more general case of nonlocality [38] . On the 
contrary, it follows from the structure of the polarization 
operator in a magnetic field [Eq. (10)]. 



B. The double trace 



The double trace of Eq. (167) turns out to be 



(kfY 



ftiy 



-6£5 + 2J(£ 0(S -% 5 ) 



(168) 



This is, in general, a nonvanishing quantity. This 
statement can be verified, for instance, by considering 
the weak- field approximation of the Euler-Heisenberg 
Lagrangian. In QED, "weak" means small as com- 
pared to the Schwingcr characteristic value B c = 
m 2 /c = 4.42 x 10 13 G, provided J > 0. In this asymp- 
totic regime the field derivatives read [39] 



5 45ttB2> 



4a 1 
^ ~ 457sF ■ 



£00 



7a 1 



(169) 



Then the double trace of Eq. (168) becomes 



(kfY 



It is worth mentioning at this point that the double trace- 
lessncss condition is customarily taken in GLVE on the 
grounds that the trace may be absorbed into the field 
renormalization [9]. In QED, the (infinite) renormaliza- 
tion has been fulfilled at the stage of the one-loop cal- 
culations that underlie the QED expressions for £ and 
for its field derivatives, the polarization tensors. The 
standard renormalization procedure of QED relates only 
to the zero-background field limit, while the background 
field-dependent part is fixed and obeys the condition re- 
flecting the correspondence principle: £5 — > as B — > 0. 
This condition is respected by Eq. (169) and establishes 
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the absence of radiative corrections to the Maxwell La- 
grangian for small and steady fields. So, the Maxwell 
Lagrangian remains untouched in this limit, which is the 
physically necessary requirement. Therefore, no renor- 
malization of the electromagnetic field additional to the 
one performed in the course of the infinite renormaliza- 
tion procedure is admitted. All the terms in Eq. (168) 
are physically important, as they serve various compo- 
nents responsible for the vacuum refraction processes 
[Eq. (40)]. We shall see in Sec. VID that this situation 
is retained even in the nonbirefringence case. 



C. Magnetoelectric coefficients 

Let us introduce, as is customarily, the matrix 
combinations (kdeY 3 = — 2 (kf) 0i0j , (khbY j = 



i e iki e j Pq , Kp 



klpq 



and (k db ) 13 



-(khe. 



ji 



kpq 



(k f ) 



Ojpq 



In terms of these quantities, the most gen- 



eral form of the quadratic Lagrangian is [10] 



2 



5 ij + (kdeT 



2 



5 ij + {khbT 3 \ bj 
(170) 



We wish to specialize this expression to the case where 
a magneticlike background (# > 0, © = 0) induces LSB. 
To this end, we compare Eq. (170) with Eq. (14). As a 
consequence, the following relations are established: 



= + {kde) 

(KDE)ij = ~ 



i J • 



(171) 



The derivation of these relations requires the use of the 
corresponding optical tensors given in Eq. (39) and is in 
agreement with Eq. (167). The remaining matrices, i.e., 
(kdbY 3 and (khe) 31 , are responsible for a magnetoelec- 
tric effect, which is the magnetic linear response to an 
applied electric field and, reciprocally, the electric linear 
response to an applied magnetic field. These matrices 
vanish identically in our framework. For the magneto- 
electric effect, and hence for the matrices (kdbY 3 and 
(khe) 31 '> to exist, it is necessary to admit [5] the nonvan- 
ishing of the pseudoscalar invariant of the external field, 
0^0, i.e., to take a general combination of an electric 
and a magnetic field as the external field. 



D. The anisotropic nonbirefringent case 

The special option of anisotropy without birefringence 
is often paid attention in GLVE; for example, in Ref. [40]. 
To establish the conditions for the absence of birefrin- 
gence in our context, we must equalize the dispersion laws 
for the two different cigenmodes of Eq. (12), = 



fa(kj_). It follows from Eq. (38) that in the infrared 
limit, one has [28] 



Mkl) = kl 



l-£a 



h{k\) = kl 1 



(172) 
(173) 



Hence, there is no birefringence, provided that the ef- 
fective Lagrangian is subject to the condition 



23 , £©e5£s5 —(1 - £y) (£ && — £55) 



(174) 



This condition is Lorentz invariant: once there is no bire- 
fringence in a special frame, there is none in any iner- 
tial frame. Note that the nonbirefringence condition is 
not the condition of coincidence of the eigenvalues >ci 
and ^3 [Eq. (38)] (which does not take place even on 
the common mass shell of the two eigenmodes [28]), nor 
the coincidence of the two dielectric permeability values 
[Eq. (40)], contrary to what one might think. The con- 
dition of Eq. (174) is not fulfilled in QED, where the 
Heisenbcrg-Eulcr Lagrangian [34-37] is taken for £. The 
only Lagrangian where the background field tensor makes 
up its single argument free of birefringence is (as inferred 
in Ref. [41]) the Born-Infeld Lagrangian [42], 7 wherein 



BI 



1 



2^ 

„2 



3/2 



1/2 



, £ 



BI 
00 



1 
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„2 
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Here, a is the dimensional parameter inherent in that 
model, and the requirement of the correspondence prin- 
ciple £#|s-_ = is obeyed. In the small external field 
domain the two dispersion curves [Eqs. (172) and (173)] 
become 



h {k\ 

h (kl 



Iff-X) 



= fci (1 - 2ff£ ee ) , (175) 
= *i(l-2ff£ w ). (176) 



Therefore, in this domain, the nonbirefringence condition 
[Eq. (174)] reduces to £®@ = £$$\$ =0 and to the disap- 
pearance of the last (Weyl-like) term in Eq. (167), the 
same as in Ref. [40]. This fact allows us to represent 
Eq. (167) in the same form as in GLVE [40]: 



'r) Kfi K Xv 



?] X ^K KU 



n Xu K Kft ) 



(177) 



However, the symmetric tensor k^ v contained in this ex- 
pression 



K^ = (KF)^ U +W 
1 



=■ (-% + 2ff£ w ) rT + % S ^ A ^ A ^,(178) 



7 Another example of the absence of birefringence is supplied by 
noncommutative electrodynamics in an external field [43] 
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where (kf)\ v is the trace of Eq. (167), is not traceless. 
Its trace, k£J = — 2£y, disappears only in the zero exter- 
nal field limit, but otherwise, in the small-field regime, 
it depends quadratically on the external field strength- 
the same as other terms in Eq. (167) or in Eq. (178)- 
and forms the whole of the isotropic part of the dielec- 
tric and magnetic permeability tensors in Eq. (39). In 
other words, it determines the isotropic part of the vac- 
uum polarization in the external magnetic field. This is 
a physically meaningful quantity that cannot be expelled 
from the Born-Infeld model, the same as the double trace 
[Eq. (168)] from QED. 

The case under consideration is parameterized by two 
quantities £5 and £gj, with the external magnetic field B 
driving LSB. There exists, however, an additional param- 
eter K tr = IS'fijj — £5, whose definition is intrinsically 
associated with the matrices that characterize the the- 
ory [see Eq. (180) in the next subsection]. We combine 
the latter expression with the convexity properties of the 
effective Lagrangian [Eq. (41)] to establish the condition 8 

Ktr > -1. (179) 

Let us finally remark that the relation above must be 
understood as a direct consequence of fulfilling the fun- 
damental unitarity and causality principles. 



E. Numerical estimates 

We now estimate to what precision the coefficients 
(kf) k ^ Xv should be measured in order that the Lorentz 
violation caused by the magnetic field of a given mag- 
nitude might be detected in the vacuum, and con- 
front it with sensitivities attained in existing experi- 
ments aimed at detecting the intrinsic Lorentz viola- 
tions within GLVE, as these are listed in Ref. [44]. 
To this end, we begin with the matrix combina- 
tions {ke+Y 3 = 5 [(kdeY 3 + (kHbY 3 ] , (Ke-Y 3 = 

\ [(kdb)' 3 ' — (khbY 3 ] - \S ij (K DE f, and {k ±f = 
\ [(k-dbY 3 i ( k heY 3 ]i which are frequently used in de- 
termining the parameter space to which the experiments 
on birefringence are sensitive. We find it convenient to 
remark that these refer to the context of GLVE. The sub- 
stitution of Eq. (171) into the latter set of matrices allows 
us to to express {k. ±Y 3 = 0: 

1 2 

(k e+ Y 3 = -fyS* + \ [£«« - £ SS ] B l B 3 < 1(T 32 , 

(K e - Y 3 = Kt T S ij + i + £55] B l B ] < 10- 17 . 

(180) 



Inequalities [Eq. (41)] are fulfilled in the Born-Infeld model 



Here the inequalities indicate the experimental sensitivity 
related to the measurement of the corresponding coeffi- 
cient combinations. 

In considering a possible magnetic field-like Lorentz 
symmetry violation, we must restrict ourselves to very 
small magnetic fields. With QED in mind, we should 
then refer to Eq. (169) in this case. Note that 
^ = 1.03 x 10~ 4 . The accuracy of 10" 32 would be 
enough to detect a Lorentz violation produced by the 
magnetic field on the Earth, -B oart h = 0.3 — 0.6 G; the ac- 
curacy of 10~ 14 would allow one to fix a Lorentz violation 
equivalent to the presence of a magnetic field on the order 
of 10 9 G. This is an unearthly large value, on the pulsar 
scale. However, an increase in experimental accuracy- 
say, to 10~ 21 would give us the possibility of detecting 
the electromagnetic wave refraction due to an external 
magnetic field already at the laboratory value of 10 5 G. 
That would be great, because up to now no effect of bire- 
fringence of QED has been seen in the vacuum in a direct 
experiment [45, 46]. Unfortunately, as long as any mate- 
rial strongly enhances the effect of the magnetic field as 
compared to the vacuum, the above considerations may 
only relate to experimental devices that do not exploit 
matter as a medium for electromagnetic wave propaga- 
tion. Apart from such devices, electron-positron pair cre- 
ation by a single photon, and the photon splitting and 
merging are the well-recognized, efficient processes that 
are efficient in pulsar magnetospheres at magnetic field 
strengths above 10 12 G. The exclusion of one-photon pah- 
creation at the accuracy level of 10~ 20 from an ultrahigh- 
energy cosmic ray event in Ref. [40], if viewed within 
within QED, only implies that there is no magnetic field 
larger than 10 6 G in the space region where that event 
occurred anyway, which is not unexpected. 

The general conclusion of this subsection is that the 
prospects of detecting Lorentz violations in the vacuum 
by perfecting the existing experimental means might be 
based only on the belief that these violations are for some 
reason much larger than the ones induced via QED by 
magnetic fields (presumably present in the Galactic back- 
ground). 



VII. SUMMARY AND OUTLOOK 

The main effort in this work has been to perform an 
analysis of Lorentz symmetry breaking by an external 
field in nonlinear electrodynamics, the gauge sector of 
QED included. We defined how transformations from 
the residual symmetry space-time subgroup, left after the 
external time- and space-independent magnetic field had 
been imposed, act on coordinates and other vector enti- 
ties via the external field tensor. For small and steady 
electromagnetic excitations over the magnetic field back- 
ground, we have developed the Hamiltonian formalism 
to serve the linear electrodynamics of the equivalent 
anisotropic medium, for which purpose a quadratic La- 
grangian of these excitations is written with the help 
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of the polarization operator in the external magnetic 
field. The electric and magnetic permeability tensors in 
this linear electrodynamics are shown to be those of an 
equivalent uniaxial medium in any special Lorentz frame, 
where the external field is purely magnetic or purely elec- 
tric. Their principal values are expressed in the paper in 
terms of the field derivatives of the effective Lagrangian. 
The fields and inductions are given the sense of canoni- 
cal variables, the necessary primary and secondary con- 
straints are determined, and the Dirac brackets are de- 
fined on constrained physical phase space in accordance 
with the U(l) gauge invariance of the theory. The con- 
served Noether currents, corresponding to the residual 
symmetry transformations, as well as the nonconserved 
Noether currents corresponding to the Lorentz transfor- 
mations, the symmetry under which is violated by the 
external field, arc defined on the physical phase subspace 
of the problem. Among the former, and of some im- 
portance, is a nonsymmetric but gauge-invariant energy- 
momentum tensor, used to form the Poynting vector and 
the momentum density, which are not the same quan- 
tities due to the antisymmetric part of the tensor. Wc 
have calculated the Dirac bracket commutation relations 
between all the generators of infinitesimal space-time ro- 
tations and translations to see that the SO(3) algebra of 
the photon angular momentum remains intact, despite 
the violation of the rotation symmetry by the external 
field, whereas the Poincare algebra is distorted. We de- 
rived the evolution equation for the photon angular mo- 
mentum, which is governed by the photon magnetic mo- 
ment depending on the antisymmetric part of the energy- 
momentum tensor. The polarization effects entering the 
equation of motion of the photon angular momentum 
are closely associated with the existence of an optical 
torque. This is a phenomenon inherent to conventional 
electrodynamics in anisotropic media, which manifests 
the breakdown of the rotational invariance. We argued 
that a small-amplitude electromagnetic wave propagat- 
ing in a strong magnetic field behaves as a quasiparticlc 
carrying a gauge-invariant magnetic moment orthogonal 
to the wave- vector. The corresponding analysis of the 
equation of motion for the angular momentum of light 
in a weak magnetic field was not developed here. This 
limiting case seems to be very convenient for probing the 
nonlinear behavior of the quantum vacuum. The latter 
could be achieved by transferring the angular momentum 
from the small waves to a microscopic absorptive object 
(e.g., tweezers). However, a detailed analysis of this issue 
will be given in a forthcoming publication. 

It would be also a challenge to find a closure to the dis- 
torted Poincare algebra for the present case-or, perhaps, 
other, simpler, cases-of the Lorentz-symmetry violations. 

In the last Sec. VI, we confronted the structure of the 
polarization operator in the magnetic field with prescrip- 
tions of a general Lorentz-violating electrodynamics, and 
discussed the common features and differences. 
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Appendix A: TWO IMPORTANT BRACKETS 

1. The case of {/, ,<3* } f 

The aim of this appendix is to compute the Dirac 
bracket {^,& >0 } x - In order to obtain the latter, we 
divide it into three terms: 

i ago 



i a ago 



i x ago 



(Al) 



J z = -Jd 3 yd 3 z{J i ,<p a (y)}C-i(y,z) 

An explicit derivation of the expressions above requires 
us to know the functional derivatives associated with the 
photon angular momentum [Eq. (113)]. For further con- 
venience, we write the latter as 

J = J d 3 x[-(x -b)ir + (x ■ir)b+ (x x a)V • tt](A2) 

In correspondence we obtain 
S 



Sir 1 



-(x-b)5 lj + Vx l -V l (xxa) j , 



(A3) 



5a l 



- e jlk V k (x • tt) + (x x V) V - e jlk x k V ■ n. 



The explicit structure of can be found in Eq. (65). 
Because of Eq. (93) wc will ignore the contribution pro- 
portional to ao- Having this in mind, the respective 
derivatives turn out to be 



5 0>O 



*H ^ = (Vx/>)<. (A4) 

Substituting Eqs. (A3) and (A4) into J^i 2 , one finds 

J>{ = J d 3 x [-e ijk e j V k (x • tt) - e j {x x V)% 1 

+ e ijk x k e j V ■ tt] , (A5) 

Jl = J d 3 x [6V'(V x h) 1 - (V x h)\x ■ b) 

-(Vx/ijVfixa)']. (A6) 
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The last term in Eq. (A6) vanishes identically provided 
an integration by parts. 

We expand the derivative in the first integrand of J?± 
and use the vectorial identity 

u x [(as x V) x o] = u l (x x V)o' — u ■ (as x V)o 
= -e ljk u j x k {V ■ o) + e ijk u j x m V k o m . (A7) 

As a consequence. 

Ji = J d 3 x{-(exTrf-e l (xx V)V}. (A8) 

Let us turn our attention to ■ Note that an integration 
by parts in the second term of the latter leads to 



2. The case of {JKT, &>°} <t 

Let us consider the Dirac bracket involved in Eq. (128). 
To determine the latter, we express the photon boost in 
Eq. (114) in the following form: 



,JC = Jt a? -+- JC qoo , 



JT^ = -x &>, JTqoo = / d 3 x(xQ 00 ^ 



(A15) 



J d 3 x [&V(V x h) j + (hx bf + e ijk h k x l V j b l ] . 

(A9) 



Adding to this expression a vanishing term 
- J d 3 xe tjk h j x k V ■ b and using Eq. (A7), we find 

= J d 3 x[(hxby - h\x xvyb 1 ] . (aio) 

Now, we focus ourselves on J^. Expanding the sum 
over a and /3 and taking into account Eq. (84) , it reduces 
to 



where S 1 is the spatial translation generator and 00 is 
the energy density see Eq. (65). As before, we neglect 
the contributions involving the Lagrangian multiplier ao- 
Considering Eq. (A15), one finds 

{jr,^* }^ = {.jr^,^ }^ + {Jtr£oo,&°} t (Ai6) 

The first bracket on the right-hand side can be com- 
puted by applying Eq. (135). According to this equation, 
{^\^ }^ vanishes identically. Thus 



{JTJ,,^ }^ = -x°{ 



S r7}>i raO 



}„ = 0. (A17) 



In order to analyze the remaining terms on the right- 
hand side of Eq. (A16), it is rather convenient to have at 
our disposal the following set of derivatives: 



5lT l 



-x e , 



00 



x\V x h) 1 



(A18) 



'3 = 



J d 3 yd 3 zVl {J\MV)} ^i S{3) (y- z )V Z m e mn 
x {a„(z),^°}~ J d 3 yd 3 zVy m e mn {j\a n (y)} 

(All) 



With the above expressions in mind, one can undertake 
the calculation of the second bracket in Eq. (A16). We 
write the latter as 



x ^-5^(y-z)Vl{n k { Z y.^} 



{^oo,^ h = {j^oo,^°} + ^ 



(A19) 



With the help of Eqs. (119) and (96), it is simple to 
obtain the following identities: 

Vl{/>fcW}=0, V z fe {^ fe (z),^°}=0. (A12) 

In correspondence, vanishes identically as well. We 
then substitute Eqs. (A8) and (AIO) into Eq. (Al). As 
a matter of fact, 

{/.^°}, = / d 3 x{{x xVyO 0Q +ir xe + hxb} , 

(A13) 

where the identity 

(x x V) 6 00 = -e'(as x V)tt' + h\x x V)6'(A14) 

has been taken into account. Note that an integra- 
tion by parts carries out a change of sign in the in- 
tegral / d 3 x (x x V) 9°° = - / d 3 x (x x V) e 00 . There- 
fore, the latter vanishes identically and can be ignored 
in Eq. (A13). This operation leaves us with the terms 
which appear in Eq. (126). 



with 



•W 



d 3 yd 3 z {j4oo,<p a (v)} Cj(y,z) {^(z), . 

(A20) 

From Eqs. (A18) and (A4), it is straightforward to get 

{JT,^ } = Jd 3 x{exh)\ (A21) 
Withthe help of Eq. (84) , we express W as 
W = Jd 3 yd 3 zVi{,J^ 00 ,n k (y)}^^6^(y~z)S7 z ml 
x {a„(z), ,^ } - / d 3 yd 3 zV y m e mn {j^ 

(y)} 



(A22) 



Observe that the last integral vanishes identically be- 
cause V% {-Kkiz),^ } = [see Eq. (A12)]. The first 
integral vanishes identically as well because 

Vn^eoo, 7 r fc ( J/ )}=V? d ^ 



Sa k (y) 

x l V ■ (V x h) = 0. (A23) 
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We then conclude that W — 0, and the Dirac bracket Furthermore, the substitution of Eqs. (A3) and (A18) 
is just as it appears in Eq. (129). into X 2 yields 



Appendix B: DERIVING THE ALTERATIONS TO 
THE LORENTZ ALGEBRAIC RELATIONS 



1. Dirac bracket between X and ^ 

In order to determine the Dirac bracket {J(f, we 
consider Eq. (A15). This allows us to express 

S*}. = {4, «/''}, + • (bi) 

The first bracket on the right-hand side can be computed 
by applying Eq. (141). Indeed, according to this equa- 
tion, = e l]k ^ k provided that the fields van- 
ish at infinity. Therefore, 



{-Xkf j }* = z ijk [-x 



askl 



(B2) 



Analogously to Eq. (Al), we split the second Dirac 
bracket in Eq. (Bl) into three different contributions: 



{j£goo !t / J '},=Xi-X 2 +X 3 , 



(B3) 



with 
Xi = / d 3 x 



X 2 = J d 3 x {-zWv^as • x) - x i e\x x V)V 

+ xV'Va^V -x)} . (B8) 

Nevertheless, the desirable expression of Z 2 is obtained 
by expanding the derivative present in the first term and 
using the identity given in Eq. (A7). With these details 
in mind, we find 

X 2 = J d 3 x{-x l [e l (x x V) J Y + (e x x) J ] } . (B9) 

On the other hand, once X3 is expanded over a and /3, 
one obtains 

X3 = J d 3 yd 3 zVt{j^ 00l n k (y)}^^S^(y-z)V* m s mn 
x {a n {z),J J }- [ d 3 yd 3 zyy n e mn {je^ 00 ,a n (y)} 

(BIO) 



Sai 5tt 1 ' 



Sni da 1 ' 



(B4) 



X3 = - J d 3 yd 3 z {JT, <p a (y)} C~l (y, z) Wp(z),J) 

The explicit insertion of the derivatives involved in X\ 
allows us to write 

X x = J d 3 x [e ijm h m {x ■ b) - e ilm h m Vx l -i'(Vx h) j 

x (x ■ b) + x l (V x h) l b>x l - x*(y x h) l X7 l (x x ay 
+e Um h m V{xx a y] (B5) 

Thanks to the vectorial identity given in Eq. (A7), we 
are able to express the integral of the first four terms in 
this equation as 

- J d 3 x [x l (h\x x Vyb 1 + (h x b) 1 )] (B6) 

Note that the fifth and sixth terms in Eq. (B5) cancel 
each other, since an integration by parts leads to 

J d 3 xx l {V x h) l V\x x a) 1 = J d 3 xe ilm h m Vi(x x a y . 

Once this fact is taken into account and Eq. (B6) is in- 
serted into Eq. (B5), its right-hand side acquires the fol- 
lowing structure: 

Xi = J d 3 x{-x l [h\x x V) j b l + {hx by] } . (B7) 



where the relevant elements of [Eq. (84)] were in- 
serted. Thanks to Eq. (119), both integrals in I3 vanish 
identically, and one ends up with X 3 = 0. Hence, by 
substituting Eqs. (B7) and (B9) into Eq. (B3), we get 

- J d 3 x {x l [(ir x e) 3 + (hx b) 3 }} (Bll) 

where Eq. (A15) is taken into account. We then use the 
definition of t^ v [Eq. (132)] to write the second line of 
this equation as it stands in Eq. (142). 



2. Dirac bracket between J? % and ^ J 

We start off by expressing the Dirac bracket between 
f % and ^ J in terms of two elements 



» 03 



} sr =h+h 



h = / d 3 x 



Sa l Sir 1 Sir 1 5a l 



h = -J d 3 yd 3 z {J,<p a (y)} Cj(y,z) . 

(B12) 
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An explict substitution of Eq. (A3) into I\ allows us to 
write 

h = Jd 3 x{2e ijk [-(x-6)V fe (x-7r) + (V-7iV(x-6)] 

+ e Um b j x l V m {x ■ tt) - ^ lm b l x l \J m {x • tt) + e llm x m 

x (V •tt)V / (x x a) 1 - t llm x m {V -tt)V 1 (x x a) 1 

- (x ■ b)(x x V)V - {x ■ b)(x x V)V - e l(m V m 
x (x-7r)V z (x x a) J +e j7m V m (x-7r)V / (x x a) 1 

- V z (x x a) j {x x V)V + V'(x x a) l (x x V)V} 

(B13) 

where the terms that vanish due to the antisymmetric 
property of e* J have been omitted. Note that the last 
four terms of this expression vanish as well provided an 
integration by parts. The remaining integrant can be 
written as 

h = e' y ' fe J d 3 x{2 [-(as • 6)V fe (x • tt) + (V • n)x k (x ■ b)] 

- (x ■ tt) [(x x V) x &f - (V • tt) [(x x V) x (x x a)] fc 
+ (x -6) [(x x V) x Trfj . (B14) 

We then use the identity (x x V) x A = — xV • A + x l VA l . 

As a consequence, 

h = e' y ' fc J d 3 x{-2(x-6)V fe (x-7r) + (V -n)x k 

x (x • b) - (x • 7r)2;'V fe fe' - (V • ir)x k (x ■ b) 
- (V •*)(* • 6)x J V fe (x x a) + (x • 6)x i VV} . 

(B15) 

Expanding the derivative of the first integrand and in- 
tegrating by parts the third and fifth terms, we end up 
with 

h = e t3k J d 3 x{-(x-b)n k + (x-w)b k 

+ (V -tt)(x x a) k } , (B16) 

where Eq. (A2) was inserted in order to obtain the second 
line. 

Now, it is rather clear that I2 reduces to 

h = Jd 3 yd 3 zS7t{^\7r k (y)}^^6^(y-z)\7 z m e mn 

x {a n {z), J? j } - J d 3 yd 3 zV y m e mn {/\a n {y)} 

x -±-5^{y-z)Vl{Mz),/ 3 }, (B17) 

where Eq. (84) is used. Since V" { a f i ,'^k{y)} vanishes 
identically [see Eq. (119)], h does not contribute to the 
Dirac bracket between J? % and Equipped with this 
result and substituting Eq. (B16) into Eq. (B12), we end 
up with the bracket written in Eq. (142). 



3. Dirac bracket between Jfc 1 and Jff 3 

Let us conclude the derivation of the modified Lorentz 
algebra by obtaining the Dirac bracket . A 

straightforward substitution of Eq. (A15) leads us to 
write this as 

+ J£qoo , + |^^qoo , J^oo j(B.18) 

The first three brackets are easy to compute. For in- 
stance, 

{jfj^JTj,} = x 02 {^\^'}^ = (B19) 

where Eq. (137) is used. Also, by considering Eq. (95) 
and using integration by parts, we find that 



= x ^V'. 



(B20) 



As a consequence, 

I , J^qoo | +|^eoo,^^| =0. (-B21) 

All that remains is to compute the last bracket on the 
right-hand side of Eq. (B18); i.e., the Dirac bracket be- 
tween J^qqq and J^qoq. To compute this, we write 

|=/^qoo , ^qoo I — Si + E2, 



Si = / d A x 



S 2 = - J d 3 yd 3 z{JXT 0OO , <p a (y)} C a ^(y,z) {<pp(z), Jf e oo} 

Inserting Eq. (A18) into Si, one finds 

Si = e ijk J d 3 x { x x (e x h)] k ) . (B22) 

This bracket can be written in a more appropriate form 
by using the relation 

(e x h) 1 + (tt x by - a ; V • tt = r l °, (B23) 

where a term proportional to ao has been ignored. 
Indeed, the substitution of the expression above into 
Eq. (B22) allows us to get 

Si = -e ijk / k - f d 3 x (xV° - x j t m ) . 



(B24) 

Thanks to Eq. (84), S 2 reduces to 

S 2 = J d 3 yd 3 zVl{,je\^ k (y)}^-^5^{y~z)V z m e n 

x {a n (z),je j }- J d 3 yd 3 zV v m e mn {,X\a n {y)} 
1 



6^(y-z)Wt{-K k (z),^} 



(B25) 



2G 



But, since V| {j^ 1 , 7Tfc(j/)} vanishes identically [see 
Eq. (120)], we can assert that E 2 does not contribute 
to the Dirac bracket of Jtf % and J^K Due to this fact, 



c ijk 



d 3 x (xV° 



xV°] 



(B26) 



Appendix C: THE MOMENTUM AND PHASE 
VELOCITY OF THE EIGENWAVES 

We find it convenient to determine a connection be- 
tween the translational generator associated with each 
degree of freedom and its respective phase velocity. In 
order to do this, we express the electric field of each eigen- 
mode as it is given in Eq. (159). Likewise, the magnetic 
field of each eigenwave turns out to be 



& (A) (*0 r , i 
—, , cos^ax — k ■ x\ . 



\bW(k)\ 



(CI) 



As in Sec. V, T^ X) and oj\(k) are the amplitude and 
frequency of mode A, respectively. Besides, whatever 
the nature of the external field, the unit vectors ~ 
e( A )(fc)/|e( A )(fe)| in Eq. (159) and ~ b^ (k)/\b^ {k)\ in 
Eq. (CI) must be understood as the respective electric 
and magnetic polarizations. 9 We remark that the plane 
wave decomposition for the induction vectors d and h fol- 
lows from Eqs. (159), (CI), and (15). With these details 
in mind, the Maxwell equations (72) and (96) and (97) 
read 

Jfe • d^ (x, x°) = 0, Jb x e ( A > (x, x°) = w A 6 (A) (x, x°), 

k-b^(x,x°) = 0, A; x h^(x,x°) = -uj x d^(x,x°). 

(C2) 

Thanks to the Faraday equation, the momentum associ- 
ated with each propagation mode 



d 3 x S» x b^ 



(C3) 



can be written as 



d 3 x 



d 3 x 



u>x(k) 

1 

u\(k) 



dF> x (jfe x e< A >) 



(C4) 



where u\ = uj\/\k\ is the phase velocity and n = k/\k\ 
denotes the wave vector. 



9 When the external background is a magneticlikc field tensor, 
i.e., J? > and = 0, the behavior of e (x ~>(k) and &( A )(fe) can be 



found below Eq. (13). 



Now, the Ampere law allows us to obtain the following 
relation h w = u x (n x d (A) ) + (n ■ ft (A) ) n. Multiplying 
the latter by 6, we end up with 



h W. b W 



<t A 



nxd^l -b^ 



u x n 



(A) 



d (A). e (A)_ 



x &w 

(C5) 



Its substitution into the energy [Eq. (58)] yields ^°( A ) = 
/ d 3 x [d^ ■ e^ A )] . We use this identity to express 
Eq. (C4) in the following form: 



30(A) 



u\(k) 



n. 



(C6) 



Thus, the translation generator associated with each II^,, 
cigenmode turns out to be parallel to the wave vector. 
Observe, in addition, that Eq. (C6) allows to write the 
phase velocity as u A = /\& w \n. 



Appendix D: THE POYNTING VECTOR AND 
GROUP VELOCITY OF THE EIGENWAVES 



Let us turn our attention to the Poynting vector given 
in Eq. (54). In order to simplify our exposition, we will 
confine ourselves to the case in which the external back- 
ground is a magneticlikc field > 0, = 0). The re- 
sults, however, are easily extensible to the case of an elec- 
triclike vector (^ < 0, & = 0). To establish a comparison 
with the previously discussed translation generator, it is 
rather convenient to work with the spatial integral of 
Eq. (54): 



d 3 x Q l0 (x,x°). 



(Dl) 



It also advantageous to express the Poynting vector in 
terms of it and 6. This is carried out by inserting the 
relations e = -ir/e± + ^-(n ■ B)B and h = e±b - 

£35 (b ■ B)B into Eq. (54). Considering these details we 
obtain 

<P = + xB (D2) 
where Eq. (61) has been used, and 

jT = Jd 3 x^- (b-B)v-^ (tt.B)&}. (D3) 

We remark that the last term in Eq. (D2) does not point 
in the same direction of k. Of course, each mode has a 
spatial integral of the Poynting vector given by !?( A ) = 
x B, with 



JY^ = 0, 

^ = - 



^25 <S 



d 3 X |7T 



r(2) 



jrm = M J d 3 x{b^-B)^ 3 \ 



(D4) 

B)b (2 \ (D5) 
(D6) 
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To derive these expressions, we have considered Eq. (D3) 
and the equations below Eq. (13). One must note that, 
for the second and third propagation modes, the following 
relations hold: (tt^ ■ B)b™ = —B x (ir^ x b^) and 
(&( 3 ) -B)n^ = Bx (tt^ x6^), respectivelly. This allows 
us to express Eqs. (D4)-(D6) as 



.***B x ^(2) 



(D7) 



with ^( A ) given in Eq. (C6). As a consequence, the 
spatial integral of the Poynting vector associated with 
each eigenmode reads 



(i) = ^a(i) £>(2) = ^a(2) _ 2 S-gg« ^(2) 



11 

rp(3) = 0t(3) _ ggggg ^(3) 



(D8) 



where the energy conservation (d&° /dx° = 0) is taken 
into account. Making use of the continuity equation 
[Eq. (51)] and integrating by parts, one obtains 



fP 1 /xB 



(Dll) 



where Eq. (D2) has been considered as well, and u = 
u(k)/\k\ denotes the phase velocity of the small electro- 
magnetic wave. Obviously, the velocity of energy trans- 
port associated with each eigenwave follows from this ex- 
pression and Eqs. (D7) and (D8). In this context, 



1 

Vcml = — n, 

Ul 



1 2#£ e(S 
v C m2 = — n nj 

U2 £\\U 2 J_ 



1 



-n 



£±u 3 ± 



(D12) 



In accordance with the expression above, we can conclude 
that as long as k± ^ 0, the direction of the energy prop- 
agation in each physical mode differs from its respective 
momentum. 

To proceed in our analysis, we consider the center-of- 
mass energy associated with the electromagnetic wave. 
The latter can be defined by 



1 



d xxQ 



(D9) 



where 9 00 and s0 are given in Eqs. (52) and (58), re- 
spectively. Note that the derivative with respect to time 
allows us to define the velocity of energy transport 



dxr 



1 



dx° 



d xx 



d® 



on 



dx° 



(D10) 



With U A _L = ^W/I^j^l. 

Now, we consider the dispersion equation [Eq. (11)] 
with the infrared approximation of the vacuum polariza- 
tion tensor given in Eq. (38). The corresponding solu- 
tions are given by 



<jj 2 



u 3 



(D13) 



With these expressions in mind, it is a straightforward 
calculation to show that the group velocity V\ = duj\ j dk 
of each eigenwave coincides with the respective velocity 
of energy transport [Eq. (D12)]. 
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